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Abstract 

The Hilbert spaces of the edge excitations of several "paired" fractional quan- 
tum Hall states, namely the Pfaffian, Haldane-Rezayi and 331 states, are con- 
structed and the states at each angular momentum level are enumerated. The 
method is based on finding all the zero energy states for those Hamiltonians 
for which each of these known ground states is the exact, unique, zero-energy 
eigenstate of lowest angular momentum in the disk geometry. For each state, 
we find that, in addition to the usual bosonic charge-fluctuation excitations, 
there are fermionic edge excitations. The wavefunctions for each case have 
a similar form, related to Slater determinants, and the edge states satisfy a 
"projection rule", that the parity of the number of fermions added to the 
edge equals the parity of the charge added. The edge states can be built 
out of quantum fields that describe the fermions, in addition to the usual 
scalar bosons (or Luttinger liquids) that describe the charge fluctuations. 
The fermionic fields in the Pfaffian and 331 cases are a non-interacting Ma- 
jorana (i.e., real Dirac) and Dirac field, respectively. For the Haldane-Rezayi 
state, the field is an anticommuting scalar. For this system we exhibit a 
chiral Lagrangian that has manifest SU(2) symmetry but breaks Lorentz in- 
variance because of the breakdown of the spin statistics connection implied 
by the scalar nature of the field and the positive definite norm on the Hilbert 
space. Finally we consider systems on a cylinder where the fluid has two edges 
and construct the sectors of zero energy states, discuss the projection rules 
for combining states at the two edges, and calculate the partition function 
for each edge excitation system at finite temperature in the thermodynamic 
limit. The corresponding theory for the hierarchy and its generalizations is 
also given. It is pointed out that the conformal field theories for the edge 
states are examples of orbifold constructions. Two appendices contain tech- 
nical details. 
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I. INTRODUCTION 



The theory of the excitations at the edge of an incompressible fractional quantum Hall 
state UU has undergone extensive development since its beginnings a few years ago 
In the integer quantum Hall effect, that is, when the bulk fluid fills an integer number of 
Landau levels, the edge excitations are essentially single electrons occupying single-particle 
edge states |J , that propagate in one direction along the edge and correspond to the classical 
skipping orbits. There is one "channel" of such edge states for each filled Landau level; 
each channel can be considered as a more or less non-interacting, one-dimensional, uni- 
directional Fermi sea. In the fractional effect, the edge excitations, like the bulk states, are 
highly correlated and cannot be described by single-electron states. The basic variables are 
density fluctuations which propagate in one direction along the edge. The quantum field 
theory which describes these is a chiral Luttinger liquid |2|||. In the simplest case, that of 
the Laughlin states at filling factors 1/q, q odd, this density mode is the only low-energy 
excitation at the edge. In the special case of the integer effect at filling factor 1, this is 
equivalent, via bosonization, to the single-electron, Fermi-sea description ||. The theory 
was soon extended []2]|| to the edge excitations of the hierarchy theory [[[], which yields an 
incompressible ground state in the bulk for all rational filling factors p/q with q odd. In 
general, some chiral conformal field theory, that generalizes the chiral Luttinger liquid, is 
expected to describe the low-energy, long-wavelength excitations. It was predicted in Ref. |7j 
that this theory would in general be the same as the conformal field theory whose correlation 
functions reproduce the bulk wavef unctions, and which describes such universal properties 
of the bulk states as the statistics of their fractionally- charged excitations, as it does for 
the hierarchy states. This deep connection implies that the properties of the edge states are 
not only of interest in their own right, but they can also be used to probe the properties of 
the underlying bulk state. Effects in tunneling into or between edge states have been the 
subject of various works [§,§. 

In this paper we wish to extend the theory of edge states to cover some other interesting 
states that have been proposed and which do not fit into the hierarchy scheme. In par- 
ticular, there are (i) the Haldane-Rezayi (HR) [[TtJ state, proposed to explain the plateau 
observed at v = 5/2 in terms of a half filling of the first excited Landau level, in which 
the electrons have no net polarization; (ii) the Moore- Read (Pfaflian) state M, again for 
a half-filled Landau level but this time with spin-polarized or spinless electrons. They are 
ground states of electron systems with special short range interactions, described later in 
this paper. The nature of the observed 5/2 state remains controversial (other suggestions 
include an alternative spin-singlet state ||11|| , which we believe to be a spin-singlet general- 
ized hierarchy state), while the Pfaffian has been proposed as an explanation for a v — 1/2 



plateau in double layer systems Hl2| , though theoretical calculations do not support this 



suggestion. Instead, they suggest |TB[ that the ground state there is (iii) the so-called 331 
state. Since this is constructed as a two-component generalization of the Laughlin states, 
it is part of the generalized hierarchy [TPl, but since it can also be interpreted, like the HR 



and Pfaffian states, as a "paired" state, it will be natural to include it here. 

The terminology "paired" state must be interpreted carefully. It is an old idea that 
Laughlin's states might be generalized if the electrons are first grouped into clusters of 
m particles (such as pairs, m — 2) and the resulting objects of charge m then form a 
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Laughlin state. While this may be possible, and was apparently part of the idea of HR 
H| , it is not quite what we have in mind. In terms of the now-popular composite fermions, 



which consist of an electron plus an even number q of attached vortices, and which at 



filling factor 1/q behave as particles in zero net field |L5j, the paired states are obtained by 



forming a BCS-type paired ground-state wavefunction, rather than a Fermi sea (the real- 
space wavefunctions are given in §11). This was first pointed out in Ref. @, and has been 
discussed more recently in Refs. [IE|jr7|| . We note that it is not clear that these two procedures 
lead to equivalent states. In particular, the first idea seems to lead to a prediction of abelian 
statistics for fractionally charged quasiparticles |]12|| , while the second has been connected 
with nonabelian statistics 0. The comparison of the two procedures hinges on the question 
whether the two operations, of grouping particles into pairs, and of attaching an appropriate 
number of vortices, commute. In any case, they do both lead to the result that quasiparticles 
have charges in multiples of l/2q rather than 1/q |7]], and, in the latter procedure, to the 
existence of BCS-type (composite) fermionic excitations obtained by breaking pairs; these 
are expected to have a gap in their spectrum. In this paper, we will see that the gap for 
the latter excitations goes to zero at the edge, and the fermions appear as gapless edge 
excitations, in addition to the usual bosonic charge fluctuations. The fermions can be 
described by quantum field theories which are related to relativistic conformal field theories 
(CFTs); in the cases of the Pfaffian and 331 states, these are the chiral versions of familiar 
Majorana and Dirac fermions, respectively. 

The goal of this paper is to understand the structure of the Hilbert spaces, and the field 
theories, of the edge excitations of these paired states. Only closed systems are considered, 
with the fluid in the form of either a droplet with one edge, or an annulus with two oppositely- 
moving edges. Knowledge of these field theories provides the necessary background for study 
of the tunneling and other properties of these states, which might be useful as a diagnostic 
for the nature of the bulk ground state. 

There is some previous work on the edge states considered here. Wen [ITSI has provided 
numerical evidence for decoupled Majorana fermions at the edge of the Pfaffian state. Wen, 
Wu and Hatsugai |T^] studied the edge excitations of the HR and other states, using tech- 
niques developed by Wen and Wu |20j for applying operator product expansions in CFT to 
FQHE wavefunctions along the lines of Jj]]. They obtained analytical results for the edge 
states, but were not able to show either the completeness or the linear independence of 
their states, though the dimensions of the spaces were confirmed numerically for the lowest 
excitations. 

In the remainder of this paper we will first write down (in §11), for the Pfaffian, HR 
and 331 states, wavefunctions for edge excitations that are zero energy eigenstates of the 
appropriate Hamiltonian and manifestly allow an interpretation as decoupled systems of 
fermion and the usual boson excited states at the edge. These states appear to be linearly 
independent, and it remains to check that we have obtained all the edge (or zero energy) 
states. This is proved in Appendix A. Linear independence is confirmed, up to the eighth 
level of excited states for the Pfaffian, and the sixth level for HR, by direct construction 
in Appendix B. For the Pfaffian, this reproduces and extends Wen's |1| numerical results, 



and analogous results for the HR state by Wen, Wu and Hatsugai Our method has 



the advantage over the computer calculations of being analytic and valid for any number of 
electrons. In §111, we present a 1 + 1 fermion field theory for the neutral part of the spectrum 
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of the HR state; the SU(2) symmetry of the system is explicit in this construction. The 
exponents in the singularity in the electron occupation number at the edge are predicted. 
We also present a CFT whose correlators reproduce the bulk wavefunctions in the manner of 
0. Finally, in §IV, we consider systems on a cylinder, with two edges, which gives further 
information about the structure of the systems. This information, a complete description 
of the number of states at each excitation level in the thermodynamic limit, is conveniently 
expressed as a partition function similar to the usual grand canonical one. The results are 
interpreted using the CFT idea of an orbifold construction. We emphasize that our basic 
results are derived without use of CFT, which is needed only in §111 and at the end of §IV, 
where the results are compared with CFTs. 



II. EDGE STATES OF A DISK 
A. Edge states of the Laughlin state 

In this section we will first review what is known about the edge states of the Laughlin 
states, emphasizing points that will help us in studying other FQHE systems. We then turn 
to results for the Pfaffian, HR and 331 states. 

Considering first the interior of a system, i.e. in the thermodynamic limit where the edge 
disappears to infinity, or in a system filling a finite but closed geometry such as the sphere, 
a FQH system at a given rational filling fraction possesses by definition a unique ground 
state (except for some global phenomena in the case of surfaces of nontrivial topology) and 
a gap for all excited states of the same or higher density than this ground state. In a number 
of interesting cases (with the particles confined to a fixed Landau level, usually the lowest) 
a model short-range repulsive interaction can be found for which this ground state has a 
wavefunction that is known exactly, and is a zero energy eigenstate. Therefore, turning to a 
QH droplet in a plane, with these interactions there is always a possibility of excitations of 
zero energy which are "inflations" of the densest zero energy state of the system. Hence a 
zero energy state always involves some deformation of the edge compared with the ground 
state. 

The simplest examples of these systems are modeled in terms of Haldane's pseudopoten- 
tials 0] by the Hamiltonian (we work in the lowest Landau level throughout) 



tf = X>£?f (2.i) 

1=0 i<j 

where Vi are positive constants (pseudopotentials) and V\ 3 is the projection operator onto the 
relative angular momentum state of angular momentum I for particles i and j. The densest 
zero energy eigenstate of H Q2.1D , that is the one with the lowest total angular momentum 
(= total degree of the polynomial part) is the Laughlin state (here in the symmetric gauge) 

m 

* L ( Zl , ...,z N ) = H(zi - ztfexpi-j £ N 2 ]. (2.2) 

i<j 
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If the Laughlin state is multiplied by any symmetric polynomial in the Zi's, it will still be a 
zero energy state [21]. We will elaborate a little on this point. 
The Laughlin quasihole operator 



N 

n(*-«o. ( 2 - 3 ) 

multiplied into produces a quasihole located at w, and many- quasihole states can be 
obtained by repeated use of this operator. All such states are clearly zero energy excitations 
for ( |2.1| ), because the wavefunction still vanishes as (z{ — Zj) q as Zj — > Zj. The operator can 
be viewed as a generating function, through the expansion 

N 

Hizi -w) = e N + --- + (-wf- 2 e 2 + (-wf- l ei + (-w) N , (2.4) 

8=1 

for the elementary symmetric polynomials, 

e n = z ii ' " ' (2-5) 

l<il<i 2 <— <in<iV 

The symmetric polynomials in the Z{ form a closed set under the operations of taking sums, 
differences and products, so they form a ring with unit element 1; allowing linear combi- 
nations with arbitrary complex coefficients, they form an algebra. All such polynomials, 
and therefore all zero energy states, are obtained as linear combinations of products of the 
e n , i.e. the e n 's form a basis for the algebra of symmetric polynomials. The e n 's can be 
obtained by integrating the quasihole operator, times a suitable factor, over all w; thus all 
zero energy states can also be obtained as linear combinations of products of integrals of 
quasihole operators acting on ^l- Single Laughlin quasihole states reconstructed from the 
e n 's via (|2.4j) are not orthogonal and may be considered as coherent states, similar to those 
of a single electron in the lowest Landau level, which can be expanded in angular momentum 
eigenstates. Likewise here, the e n are operators that increase the angular momentum of the 
electrons by n. Also, as operators, they clearly commute. 

As long as w lies well inside the disk of radius ^2q(N — 1) formed by the electrons, say 
more than one magnetic length from the edge, it represents the position of the quasihole. As 
w — > oo, the disturbance in density due to the quasihole approaches the edge and eventually, 
in the limit, becomes trivial. The leading corrections to this limit are then the terms e\, 
e2, ... so that e n with larger angular momentum represent larger distortions of the edge. 
This is borne out by the energetics, if we introduce a term AM into the Hamiltonian, where 
M is the operator representing total angular momentum about the origin. Since rotations 
are a symmetry of the unperturbed problem, this term merely splits the degeneracy of the 
different angular momentum eigenstates. Then it is clear that states e n ^i have energy 
increasing linearly with n and since gapless excitations can occur only near the edge of an 
incompressible QH state, we can again conclude that these states are edge excitations. This 
can be verified in detail for the q = 1 case which is a filled Landau level 0||]. The same 
expansion for these coherent states with more than one quasihole gives for each linearly 
independent bosonic state a corresponding edge state. By the arguments above, these must 
span the full space of edge excitations. 
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One may count the number of states at each increased angular momentum, as follows. 
Note that to describe edge states, we consider the limit N — > oo, with AM = M — M fixed 
(M is the angular momentum in the ground state, and M = qN(N — l)/2 in the Laughlin 
state). Bulk states are obtained either by fixing w's in quasihole states or by applying e n 's 
with n of order N, and then N — > oo in either case. Since AM in an edge state Yl a e n J$>L 
is AM = J2a n a, the total number of states at AM is p(AM), the number of ways AM 
can be divided ( "partitioned" ) into positive integer parts whose sum is AM. However, their 
meaning is clearest if we use a different basis for the algebra of symmetric polynomials, 
namely the sums of powers 

s„ = $>? (2.6) 

i 

(these are not all independent for finite N, but must all be used as N — > oo), which are 
one body operators, and, up to constant factors, can be viewed as the positive angular 
momentum components of the change in density at the edge from the ground state, projected 
into the space of zero energy states. Thus the edge states are built up out of repeated 
applications of density operators which do behave as boson creation operators (i.e. they 
commute) for n positive. The components of the projected density with n negative act as 
destruction operators and with a suitable normalization the algebra of independent simple 
harmonic oscillators is obtained, or equivalently, the abelian (U(l)) analogue of a Kac-Moody 
algebra. This point of view is extensively discussed in 

The partition function can be obtained from Euler's generating function which is an 
infinite series in an indeterminate x: 

oo oo 

Z(x) = l+ p(^M)x AM =U(l-x n )- 1 , (2.7) 

AM=1 n=l 

where we used the binomial expansion 

(l-x")- 1 = l + x n + x 2n + .... (2.8) 

We recognize Z(x) as the statistical-mechanical partition function of a chiral Bose field, 
that is a collection of simple harmonic oscillators of frequencies nu>, n = 1, 2, . . . , if we set 
x = exp(— w/fceT) (and h — 1). For the convenience of the reader at later points in this 
paper we include a table of the partition function p(AM) for small AM: 



AM 


1 


2 


3 


4 


5 


6 


7 


8 


dim=p(AM) 


1 


2 


3 


5 


7 


11 


15 


22 



In addition to these states, we can obtain states with a net charge added to the edge, 
either by changing the electron number, which of course can only give integral charges, or by 
adding quasiholes at the center of the disk, which allows the charge effectively added to the 
edge to be a multiple of 1/q. The wavefunctions of such a state, with a positive charge m/q 
added at the edge, are the same as the above except that a factor Hi i s included. The 
partition function for the states in each such sector is the same as for m = 0; these again 
represent density fluctuations on top of the state which now has charge m/q added. The 
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states with different electron number will, in the following, generally be found to play a role 
in the structure of the theory. Whether the states with quasiholes added at the center should 
be viewed as part of the edge theory is somewhat a matter of taste; they could alternatively 
be viewed as a part of the more general theory of bulk and edge excitations together. 



B. Edge states of the Pfaffian state 

The Pfaffian state 0, for even particle number N, is defined by the wavefunction 

y Pi (z u ■ ■ ■ , z N ) = pf ( Y[( Zi - Zj y exp[~ J2 N 2 ] ( 2 - 9 ) 

\ Z i Z jJ i<j 4 



where the Pfaffian is defined by 



! N/2 



PfMij = 2 N / 2 (N/2)\ ^ s ^ na U M ^k-iM2k) 
for an N x N antisymmetric matrix whose elements are Mj,-. It is the lowest angular 



momentum ground state of the Hamiltonian [12 



H = V ^(z.-z^iz.-Zk) (2.10) 

<ijk> 

(where the sum is over distinct triples of particles) for the case q = 1 and of similar three- 
body short-range interactions for q > 1. The filling factor is 1/q. The Pfaffian state is 
totally antisymmetric for q even, so could describe electrons, while for q odd it describes 
charged bosons in a high magnetic field. Zero energy quasihole excitations correspond to 
increasing the flux inside the area spanned by the fluid, as usual, but in this case the basic 
objects contain a half flux quantum each and must be created in pairs. A wavefunction for 
two quasiholes was proposed in 0; it is 

*(zi, . . . ,z N ;w 1 ,w 2 ) = 

1 V- nS[(2 ff (2H) - Wi)(z a{2k ) ~ W 2 ) + (Wi <-> W 2 )} 

} sgncr- 
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It is clear that for q = 1 the quasihole states are zero energy eigenstates of H, (|2.10|) ; this also 



holds for the appropriate 3 body H for q > 1. It will be seen that it is the pairing structure 
built into the ground state which allows insertion of Laughlin-like factors which act on only 
one member of each pair and hence effectively contain a half flux quantum each, unlike 
the usual Laughlin quasihole that corresponds to a full flux quantum. The same structure 
requires that quasiholes are made in pairs, since the wave function must be homogeneous. 
When quasiholes coincide, that is, when w\ = w 2 , a Laughlin quasihole is recovered. 
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The multi-quasihole states can be used to generate the edge spectrum of the Pfaffian 
state. We initially used such an approach, but now find it simpler to write down an ansatz 
which, we believe, in fact describes all the zero-energy states. We construct wavefunctions 
for N electrons (N odd or even) which we will interpret as having F fermions created at the 
edge: 

^m.,...,nF\ z U ■ ■ ■ ' z n) — 

1 Hfe = l Z a(k) 

2^ sg 110 "" 



2 (iV-F)/2 (iV _ F)/2 , ^ {MF+1) _ MF+2)) . . . _ MN)) 

xnU-^fflpl-iEN 3 ]. (2.12) 



i<j 



Here < n\ < n 2 < ■ ■ ■ < rip, are a set of distinct nonnegative integers, and N — F clearly 
must be even. The sum over permutations can be divided into a sum of terms in each of 
which the unpaired electrons (those with indices er(l) through o~(F) in the above expression) 
are antisymmetrized among themselves. Then each term contains a Slater determinant in 
these coordinates, representing fermions in wavefunctions z Uk , hence the stated conditions 
on the nfc's. The angular momentum of the states is 

F 1 

M = £ n k + -(qN(N - 1) - (JV - F)). (2.13) 
k=i 1 

Hence the angular momentum relative to the ground state, AM = M — M Q , is 

AM = f> fe + l). (2.14) 
fe=i z 

Note that the angular momentum of the ground state is calculated at the same number of 
electrons, and is M = N(q(N—l) — l)/2. Such a ground state only exists for N even, but we 
use the formula as an interpolation for N odd also, to yield ( |2.14|) . We interpret the states 
as having fermions created in orbitals of angular momentum AM = n + 1/2, n — 0, 1, . . . , 
which is exactly the description of the ground state (antiperiodic) sector for Major ana- Weyl 
fermions on a circle. We should note, however, that if we choose a certain even number 
iVo of electrons and use the ground state for iV = No as a reference ground state, then the 
odd fermion number sectors occur only when some charge has been either added or removed 
from the ground state. More precisely, odd fermion numbers can occur when N is odd, so 
an odd amount of charge has been added to the edge, and similarly for even numbers. Thus 
the parity of the fermion number is equal to the parity of the (integral) charge added. This 
seemingly trivial observation indicates that the fermion and charge edge degrees of freedom 
are not completely decoupled. This is analogous to the spin-charge separation at the edge of 
the Halperin states, where although the spin and charge form separate excitations that can 
be moved along the edge independently, there are global selection rules that relate the total 
spin and charge added at the edge, in a similar way as here |[23]. This is closely related to 



the absence of any spin-charge separation in the quasiparticles in the bulk, which can carry 
spin 1/2 only if they also carry nonzero charge. These "projection rules" will be discussed 
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more extensively later, including the hierarchy states, to one of which the Halperin state is 
isomorphic. 

In addition to these states, we can also take any one of them multiplied by a symmet- 
ric polynomial in all the z^s, which is again a zero energy eigenstate. These polynomials 
represent the ubiquitous chiral bosons associated with charge excitations and need not be 
considered further at the moment. To ensure that all these states represent linearly indepen- 
dent edge excitations, we must certainly take the limit iV — ► oo when studying each space 
of angular momentum eigenstates at AM fixed and finite. As we will see below, the states 
without symmetric polynomial factors appear to be linearly independent, but a full proof of 
this, and of independence of the symmetric polynomials, appears difficult. In Appendix A 
we prove that all zero-energy states can be written as linear combinations of the forms Q2.12|) 
times symmetric polynomials. In Appendix B we indicate how we showed, for AM up to 8, 
that all these states are linearly independent. This provides rather convincing evidence for 
our simple form ( j2.12| ). 

For completeness, we also give expressions for the edge states in the other, "twisted", 
sector, where the Majorana-Weyl field obeys periodic boundary conditions. These states 
occur when an odd number of quasiparticles are present far inside the edge. For simplicity, 
we consider a single quasihole at the center of the disk. The ground state in this sector can 
be produced by taking the two-quasihole state above, dividing by w^ 2 and letting w% — > 0, 
w 2 —>■ oo. On including unpaired electrons as above we obtain 

^ni ,...,np > Zn) — 

1 v IIfc = l Zjfo ITU (z g (F+2f-l) + Z a (F+2l)) 

2(»-W(N - F)/2\ Sgna (z CT(F+1) - z a{F+2) ) ■ ■ ■ (z a (N-i) - z a(N) ) 

xn(2i-%) ? exp[-^N 2 ]. (2.15) 

i<j 

These states have angular momentum 

F 1 

M = J2n k + -qN(N-l). (2.16) 
k=i 1 

Hence the ground state, in which F = 0, has angular momentum Mo = qN(N — l)/2 and 
the angular momentum relative to the ground state, AM = M — Mo, is 

F 

AM = J> fc . (2.17) 

k=l 

Similar remarks to those above about the F odd cases apply here. These quantum numbers 
are exactly those expected for Majorana-Weyl fermions obeying periodic boundary condi- 
tions, in which fermions can be added in orbitals with angular momentum n — 0, 1, .... In 
particular, note that for F = 1, we can take n\ = and obtain a state with zero increased 
angular momentum. This state is entirely meaningful; since it has an odd number of elec- 
trons, it is not the same as the ground state (with the added quasihole), which has an even 
number of electrons. That all states in this sector have the form shown is also proved in 
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Appendix A; proof of linear independence along the lines of Appendix B has been done at 
low AM. 

Other numbers of quasiholes, either even or odd, are obtained by taking states in the 
untwisted or twisted sectors, respectively, and multiplying in Laughlin quasihole factors 
n« z™ 1 . These states are similar to those in the two sectors above but with additional charge 
at the edge, to which it has been pushed by the quasiholes at the center. 

All these states with quasiholes at the center of the disk obey projection rules similar 
to those for the states without the quasiholes. In this CctSC, clS iV varies, the charge at the 
edge runs over values which are integers plus a fixed fraction (defined mod 1, and equal to 
a multiple of l/2q). If the charge added at the edge is defined relative to a reference ground 
state at even iV = No, for each number of quasiholes at the center, then the projection rules 
are unchanged. This avoids problems of definition in the case where q Laughlin quasiholes 
have been added at the center. A more satisfactory approach will be given for the case of 
two edges, as on the cylinder, in §IV. The states with quasiholes added at the center of the 
disk can be viewed as a special case of this. One could argue that for the disk, it is more 
natural to exclude any bulk excitations, in which case there are no twisted or quasihole 
sectors. The other sectors arise only when both edges are present, or in the presence of bulk 
quasiparticles. 

For completeness, we include a table of the dimensions of the spaces of fermion states 
for low AM in the untwisted sector: 



AM 


1 


2 


3 


4 


5 


6 


7 


8 


dim 





1 


1 


2 


2 


3 


3 


5 



The entries in the table are the ones that have been verified in Appendix B. The dimensions 
of the full space of edge excitations in the untwisted, even particle number sector are found 
by convoluting these numbers with those for the U(l) chiral boson system given earlier. Thus 
all states in this sector have been verified to be linearly independent up to AM = 8. In a 
numerical calculation in |18|] only the states for AM < 5, which in fact contain a maximum 
of two excited fermions, were found. The first state with 4 fermions excited appears at 
AM = 8 (and is included in the table). 



C. Edge states of the Haldane-Rezayi state 

The Haldane-Rezayi (HR) state PI can be written in terms of the coordinates of N/2 
up-spin electrons at z{, . . . , and N/2 down-spin electrons at 

^HR{zh - - - > Z N/2> Z lf ' ' ' 5 Z N/2) = 

E t _ J. \2 71 — ^n^-^exph-EN 2 ]- (2-18) 

creSjv/2 \ Z l Z a(l)J "\ Z N/2 Z a(N/2)) i<j 



Here q is even to describe fermionic electrons, and the filling factor is 1/q. The first factor 
is of course just a determinant. The product over z^s with no spin labels attached is over 
all particles. The fact that this describes a singlet is discussed carefully in [IU|. Strictly 
speaking, the form given is an abuse of notation. The correct way to write the functions 



10 



is as a function of N electron coordinates, numbered from 1 to N, half of which have up 
spin and half down, and the permutations are over the subset of electrons of each spin. 
On including the proper sign factors, the spatial wavefunctions can be combined with the 
spinor wavefunctions of the N electrons, and then summed over all ways of choosing which 
electrons have which spin. In this way, wavefunctions that are totally antisymmetric under 
exchange of both the space and spin labels of particles are constructed. This procedure 
is standard and has been described in the literature |T(],^|,^]; it can be used to produce 
states of definite total spin. Since the construction of such states from the functions given 
below is straightforward, if tedious, it will be omitted, and we will continue to use the 
abused notation as in (|2.18| ). In [[/]] it was pointed out that this state can be regarded as 



a BCS-type condensate of spin-singlet pairs of spin- 1/2 neutral fermions that consist of an 
electron and q vortices, from which the spin-singlet property can be more easily understood. 
The HR state is the unique zero energy state at = q(N — 1) — 2 flux of a "hollow- 
core" pseudopotential Hamiltonian that gives any two particles a nonzero energy when their 
relative angular momentum is exactly q — 1 [|K]]. 

As for the Pfaffian state, excitations of hc/2e flux are expected and by flux quantization 
they should occur in pairs. In exact analogy with the Pfaffian state, the wavefunction for 
two quasiholes is: 



N/21 
N/2, 



reS N/2 \ Z 1 Z cr(l)) "\ Z N/2 Z a(N/2)) i<j 



Due to the spin-independence of the newly inserted factors acting on each pair inside the sum 
over permutations, the state is still a spin-singlet and this suggests that the quasiholes carry 
no spin. We will see further evidence of this later. The two quasihole state is again a zero 
energy eigenstate of the "hollow core" Hamiltonian where all pseudopotentials (including 
Vo) are zero except V q ~i, i.e., for two particles to interact their relative angular momentum 
should be q — 1. To see this fact, expand the inserted factors for each pair in terms of powers 
of z\ ± zK k y Due to the symmetry between z\ and zL fe x in each factor, it is easy to see that 

t I 

z k — z^, k s must occur to an even power. Thus in the complete wavefunction, the absence 

of (z k — z^ a ^) q ~ l for any k, I, and hence the zero energy property of the ground state, is 
preserved in the quasihole states. 

It is possible to write down directly the forms of all the zero energy states of the hollow 
core Hamiltonian, in analogy with those for the Pfaffian. In the untwisted sector, in terms 
of the coordinates of iVj up electrons, iVj down electrons, the wavefunctions are linear 
combinations of (omitting symmetric polynomial prefactors) 



l T\k=i( z i(k)) nk Yli=i( z pU). 

} j sgn a sgn p- 



iiii 



( N 1 F l)- a€ s N ^ ( Z i(Fi+l) Z i(F L +l)) 2 " ' ( Z l(N^) z[ p{N y )) 2 

xnfe-^-) 9ex p[-iEN 2 ]- ( 2 - 2 °) 

i<j 
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Here iVf — Ff = JVj, — -Fj is the number of unbroken pairs, and we may assume the n^'s, 
Trie's are strictly increasing as for those in the Pfafhan edge states. These functions have a 
structure similar to the real space wavefunctions for a BCS state with some broken pairs, 
that is with BCS quasiparticles added; the latter would have a similar form for the sum 
over permutations, but the factors zj Uh would be replaced by plane waves. Here, of course, 
they represent edge states, in which the fermions do behave as if they occupied plane waves 
running along the one-dimensional edge. As written, these states do not have definite spin, 
but eigenstates of S 2 and of S z can be constructed as indicated above. Since the paired 
electrons form singlets, the spin is determined by the spin 1/2 unpaired fermions in the 
sums over a and p, which behave identically to ordinary spin 1/2 fermions. Hence the 
possible spin states are determined by adding the spins of electrons in different orbitals 
(labelled by or rrik), with the only constraint that an orbital occupied with both an up 
and a down fermion must form a singlet. 

The angular momentum of the wavefunctions given is 

M = E n k + E m k + 7:(qN(N - 1) - 2(N - F T - (2.21) 

fe=i fc=i 1 

Hence the angular momentum relative to the ground state, AM = M — Mo, is 

AM = ]TK + 1) + E( m * + 1 )- ( 2 - 22 ) 

k=l k=l 

Conformal invariance ideas suggest that this implies that the edge excitations are fermions 
of conformal weight 1, not conformal weight 1/2 as for the Pfaffian state (§IIB) and the 331 
state (§IID). This will be discussed further in §111. We note that the projection rule arising 
from our states is the same as for the Pfaffian, i.e. the parities of the fermion and charge 
numbers are the same. 

From these wavefunctions we can obtain the total numbers of the low-lying edge exci- 
tations (excluding symmetric polynomials) of the HR state in the untwisted sector at fixed 
even N: 



AM 


1 


2 


3 


4 


5 


6 


dim 





1 


4 


5 


8 


10 



Again, though the count of states based on the wavefunctions given could of course be 
continued, the table has been terminated at the largest AM where we were able to verify 
linear independence directly (see Appendix B). Note that spins larger than 1 do not occur at 
these low AM. Our numbers for spin and 1 excitations agree, when convoluted with the 
U(l) numbers, with those calculated in Wen, Wu and Hatsugai |T3fl , though our results for 
N large extend to higher AM than they can attain in a small system without encountering 
finite size effects. 

The twisted sector is again obtained by including a spin-independent quasihole factor, 

this time YliJi \ z a {F +i) + z p{F l +t))i ^ n ^ e sum on permutations. The angular momentum 
of the excited states relative to the ground state in this sector is 
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AM = + 1/2) + £(m fc + 1/2). 

k=l k=l 



(2.23) 



Finally, it is once again possible to multiply in factors IL Z T that add charge to the edge, 
which are spin independent and identical to those for the Laughlin states. 



D. Edge states of the 331 state 

The 331 state is just one of a family of two-component states, the so-called mm'n states, 
first introduced by Halperin f26|. Using notation f, J, for the two components, even though 



they need not represent spin, and bearing in mind that similar remarks to those at the 
beginning of §IIC about constructing totally antisymmetric wavefunctions apply here also, 
these states can be written: 

• • • , 4 /2 ) = 11(4 - 4) m 11(4 - 4r' 11(4 - 4) n ?M-\ E N 2 ]- (2-24) 

i<j k<l rs 4 t 

The general mm'n state is the unique lowest total-angular-momentum ground state of a 
spin-dependent pseudopotential Hamiltonian, that generalizes ( |2.1| ) to the two-component 
case, which gives positive energy to any state in which two f or j particles have relative 
angular momentum less than m or m', respectively, or in which an j and a J, particle have 
relative angular momentum less than n. 

For the case when the exponents in these states are of the form m = m' = q + 1, n = q — 1, 
q > 1, (which give filling factor v — 1/q, and the partial filling factors for f , j are both l/2q; 
for brevity, we will continue to refer to this class of states with general q as the 331 state), 
then use of the Cauchy determinant identity 

11(4 - 4) 11(4 - 4) 11(4 - 4)" 1 = det ( - T ^- r ) (2.25) 

i<j k<l rs Zj J 

allows the ground states to be written in a paired form, similar to the Pfaffian and HR 
states 0,0. This identity can be understood physically, in terms of the description of 
bulk fractional quantum Hall effect wavefunctions as conformal field theory correlators J7|, 
as expressing bosonization of correlators of a chiral Dirac (or Weyl) field (on the right hand 
side) in terms of correlators in a Coulomb gas (or exponentials of a chiral scalar Bose field) 
(on the left hand side). In terms of BCS-type pairing, this function describes p-type spin- 
triplet pairing, with each pair in the S z = state of a spin triplet flTfi|,[r7|]. 

We will extend the fermionized description immediately to include the edge excitations 
in the untwisted sector, omitting symmetric polynomial prefactors: 



i _ nri 1 (4( fc) ) Wfc nfi 1 (4 (0 ) m ' 

2^ sgnasgnp- 1 , . tt^ r 



(^t-^t)!^ r (4 ( ^ T+ i) 



pes N 



xn^-^exphjEl^i 2 ] ( 2 - 26 ) 

i<j 
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which is particularly similar to the HR case. For the angular momentum we obtain 



M = E nk + E ™* + ?(qN(N - 1) - (N - F T - Fj). (2.27) 

fe=i fc=i ^ 

Hence the angular momentum relative to the ground state, AM = M — M , is 

AM = EK + 1/2) + E(m* + 1/2)- (2-28) 
fc=i fe=i 

This is the correct behavior for the states of a chiral Dirac (or Weyl) field, where the two 
types | and j denote particles and antiparticles. This is as expected from the general 
arguments based on the form of the bulk ground state wavefunction [|7]], which, as we have 
mentioned above, includes a correlator of this same type of fields. The projection rule is 
once again the same as for the Pfaffian. 

The edge states can be reexpressed in bosonic form as 

J^tt 9+ i,9+i,9-i (2.29) 

in which (J 7 ^) are symmetric polynomials in the T (!) coordinates only, and the numbers 
iVf, iVj of | and ! particles need not be equal. 

This system also has a twisted sector obtained in a similar way as in the other examples, 
by multiplying by factors n z% or Yl z\ that represent the elementary quasiholes located at 
the center of the drop. For a single such factor, this leads to a formula for AM like ( |2.28| ) 
but in which the 1/2's in the expression are dropped. 

In the bosonic form, it is easy to see that all these states, both twisted and untwisted, 
are zero-energy eigenstates for the above-mentioned pseudopotential Hamiltonian, as an 
extension of the arguments for the one-component Laughlin states, and that they span the 
space of such states. The equivalence of the bosonic and fermionic forms of edge state 
wavefunctions involves generalizations of the Cauchy determinant identity (|2.25|) . We will 
return to the bosonized description in §IVD. 



III. FIELD THEORY OF THE HR STATE 

A. Field theory of the edge states of the HR state 

We have seen in the previous section that, apart from the charge fluctuation excitations, 
the edge excitations in the states we have studied are free fermions (for most of this section, 
we ignore the projection rules exhibited in §11; they will be reincorporated in §IV). For the 
Pfaffian and 331 states, these can clearly be described by relativistic Fermi fields in 1 + 1 
dimensions (i.e. distance along the edge and time) of scaling dimension 1/2, which we will 
denote ip f° r the Majorana field in the Pfaffian case, and if)^, ipi for the Dirac field and its 
adjoint in the 331 case. These standard field theories need not be described here. For the HR 
state, a natural candidate might have been the Dirac theory, with up and down excitations 
described by particle and antiparticle. However, in the Dirac field theory, there is no SU(2) 
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symmetry that can be generated by local expressions for the spin density and current, and 
in fact it correctly describes the edge of the 331 state. Moreover the angular momentum 
quantum numbers show that the field for the HR state does not have scaling dimension 1/2, 
but instead dimension 1. This puzzle will be addressed in this section. Another attempt at 
its resolution has been made by Wen and Wu [PD[ , which described the bulk wavef unctions, 
but did not exhibit the simple Lagrangian description shown here. 

First we write a Hamiltonian that reproduces the angular momentum eigenvalues already 
found. Introducing a velocity v s for the spin excitations when a term AM is present, the 
Hamiltonian for an edge of circumference L would be 

oo 

H = v s J2 M4t a fct + a li a H)- 

?1=1 

Here the operators a^, a\, a , with k = 27m/ L, obey the canonical anticommutation rela- 
tions {al a , cifcv} = Skk'Saa', {flfca, ctfc'cr'} = {oL) a I' CT '} = 0- Comparing with the result for 
the Pfaffian state, where there is a real (Majorana), right-moving (Weyl) fermion, (see for 
example J2[) we see that apart from the extra S z quantum number, the boundary condi- 
tion is periodic in the ground state (untwisted) sector here, while it was antiperiodic in the 
Major ana- Weyl system. 

Therefore we propose a new 1 + 1 dimensional fermion theory for the neutral part of the 
HR edge, with a doublet of complex Fermi fields ^ a (x, t), a — f or j, and a chiral Lagrangian 
density (inspired by that for chiral scalar bosons ||) of the explicitly SU(2) invariant form: 

£ = ~e™'(d t - v s d x )^ a d x ^ a > + h.c. . (3.2) 

Here £ aa ' = — e a ' a , = 1. The canonical procedure leads to the following, simpler looking 
Hamiltonian 

H = ^ J dxvsid^d^i + d^l^]) (3.3) 

together with canonical momenta = d x ^i, Hi = —d x ^fp Using periodic boundary 
conditions, and going to Fourier modes, we see that for the zero wavevector modes we 
obtain the first-class constraints that the corresponding momenta vanish: H a (k = 0) = 0. 
The constraints can be included by simply omitting the zero modes henceforth in this chiral 
theory. Quantization using canonical anticommutation relations then leads to quantized 
fields: 

^ T = ~p= ex P ~ ik ( x + v st) + E 4=^= ex P ik ( x + ( 3 - 4 ) 
k>o ykL k>0 ykL 

^j. = E 7= ex P —ik(x + v s t) — ^2 ~7== ex P + v st) (3-5) 
fc>o v kL k>0 y/kL 

and their adjoints where {al a ,a k > a >} = <W<W, {a>ka,a>k'a'} = {4<7>4v} = °> and k = 
27rn/L with n = 1, 2, . . . for a system of circumference L. It is assumed that the vacuum 
|0) obeys a ka \ti) = 0. The normal ordered version of ( |3TB| ) then yields ( |3.1[ ). 
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Because of the Fermi statistics chosen for the field, and the positive-definite norm im- 
posed, as usual, on the Hilbert space, this field theory is not Lorentz invariant, in spite of 
the gapless linear spectrum. Consequently, it is not conformally invariant either. This may 
be surprising since we have become used to the edge theories being some chiral conformal 
system, but in fact, since we started with a nonrelativistic system of electrons in a high mag- 
netic field, nothing guarantees that the edge must be Lorentz invariant, even when there is a 
linear dispersion relation for the excitations. Nonetheless, we will see that there is a closely 
related conformal field theory. 

Returning to the chiral theory, the SU(2) currents can be found using the standard 
Noether procedure: 

* = \%wkr--*' +h -* (3 ' 6) 

where we specify the Grassmann derivative 6C/5(d t ^i) to be taken from the left, r a (a = 
1,2,3) are Pauli matrices, and normal ordering is assumed. Then the total spin operators 
for the edge are S a = J dxj'^x), and 

s " = Z^(4 T a *t -°li a ki) ( 3 - 7 ) 
s + = J2 a U a H ( 3 - 8 ) 

k 

S~ = Y, a H a k1 (3-9) 
k 

which are easily seen to satisfy the SU(2) commutation relations. However, unlike other 
FQHE systems, such as the Halperin state, where the edge theory is not only conformally 
invariant but also has a Kac-Moody current algebra as a spectrum-generating algebra, here 
the currents do not form a Kac-Moody algebra, and their correlation functions contain 
logarithmic factors. 

We are ready to identify the operators describing addition or removal of electrons at the 
edge in our conjectured edge field theory. The operators cM^e -1 ^^ and d*$! ^e~ iy fi- tp represent 
the electron annihilation field operators \l/| e ; and ^ e z- in these expressions <p(x,t) is the 
usual chiral boson field representing the density fluctuations at the edge, with propagator 

(<p(x, t)(p(0, 0)) = - ln(x + v c t); (3.10) 

it is related to the electronic charge density by p = —idip/^q The exponential of the 

boson operator creates a bosonic object whenever q is even, so the scalar fermion, like the 
Majorana fermion in the Pfaffian case, field makes the whole thing into a fermion, as the 
electron should be. Note that the gradient of the scalar fermion field appears here, not the 
field itself; this reproduces the spin 1 field found earlier. Taking the charge excitations to 
propagate with velocity v c , we find that the electron propagator is 

«>l*UM)*t-(0,0)|0> oc (x + v J (x + vJ) , (3.H) 
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where the space-time separation of the two fields should be small compared with the cir- 
cumference of the disk. The total exponent is thus q + 2, in contrast to that for the Pfaffian 
and 331 state at v = 1/q which give q + l, whereas the Laughlin states give q. Consequently, 
the expectation n(k) of the occupation number of the kth single electron state, which can be 
obtained by Fourier transforming (in one dimension, along the edge) the equal-time electron 
Green's function, has a power-law singularity n(k) ~ \k — k max \ q+l for the HR state, while 
the exponent is q for the Pfaffian and 331 states, q — 1 for the Laughlin states (for the full 
Landau level, q — 1, and there is a discontinuity in n(k)). Numerical simulations have been 
performed for both the Haldane-Rezayi and Pfaffian states |Z7| , but no conclusion about the 
exponents in occupation number n(k) versus k relevant to the edge field theory is drawn in 
the published work. 

One might also expect that, just as the bulk system has pairing of the composite fermions 
0, similar algebraic BCS-type expectations should appear at the edge. Indeed, because of 
the form of the operators (j3.5|), we have, for example, 



<O|0¥ T (M)0¥ A (O,O)|O> ex (3.12) 

yX t V s Zj 

which can be viewed as a pairing function. However, this correlator omits the exponentials 
of ip needed to represent the electron operators; if these were included, the correlator would 
decay rapidly, since the fields carry the same, not opposite, charge. This could be taken 
to illustrate, for the edge theory, how pairing occurs for the composite fermions, not for 
electrons. Similar phenomena can be found in the Pfaffian and 331 states. However, while 
the correlator shown is legitimate as it stands for the scalar fermion field theory, it is not 
a legitimate correlator for the HR edge theory, because the required intermediate states, 
where only a single fermion has been added or removed from the ground state (or excited 
states, in the finite temperature case), do not obey the projection rule found in §IIC. Since 
fermions can be created from the ground state only by breaking pairs, states with an odd 
number of fermions occur only when an odd number of charges have also been added at the 
edge (assuming there are no changes in the interior of the system). Only operators that 
respect this rule can be constructed in the edge field theory. Thus, in the theory of the HR 
edge, this correlator can be constructed only for equal times, and must then be viewed as 
the expectation of a single, nonlocal, operator. 



B. Conformal field theory of the bulk HR state 

In [0 a mathematical connection between FQH wave functions and correlators in CFT 
was presented, namely the elegant wavefunctions of some important FQHE states are ac- 
tually correlators (or conformal blocks) in a chiral 2-dimensional conformal field theory. 
However, the question of what theory this would be for the HR state was left unresolved. 
The difficulty was to understand the SU(2) symmetry (in fact, the singlet nature of the 
state) in terms of CFT. We will present a solution to this problem here. The natural choice 
for the CFT in the HR case is a nonunitary theory with the Euclidian action (containing at 
this stage both right and left moving degrees of freedom for convenience): 
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j2 

S = J y^e™'d^ a d^ a , (3.13) 

and the Grassman field ty a is regarded as real. Thus is a relativistic scalar fermion 
and this model is conformally invariant, but its states do not all have positive self-overlaps 
because of the violation of the spin-statistics connection, so we say that the Virasoro rep- 
resentations are nonunitary; the central charge is c = —2. (We note that the spin-statistics 
theorem relates the statistics of the fields in a positive-definite, Lorentz-invariant field theory 
to the "spin" defined by rotations of the Euclidian two-dimensional spacetime, which in a 
conformal theory is the difference of the right- and left-moving conformal weights, and not to 
what we have been calling spin, which describes the transformation properties under SU(2) 
rotations that leave the spatial coordinates unchanged.) This system is an anticommuting 
counterpart to the system of a pair of scalar boson fields which is unitary and has c = 2. 
d^> a is a field of conformal weight 1 whose correlators 

(d* t (4)---^(i /2 )> (3.14) 

reproduce the determinant in the HR state. The action S is manifestly invariant under 
the group of real symplectic transformations Sp(2,R) (which preserves the reality property 
of the fields, and has the same complexified Lie algebra as SU(2)) and thus the correlator 
produces a singlet (since the vacuum is invariant). However the Noether currents are of 
the form ty a d^ a > — (a <-> a') which are not "good" conformal fields since their correlators 



contain logarithms, so there is no Kac-Moody symmetry. We note that Wen and Wu f20 
arrived at an equivalent description of this c = —2 CFT system in terms of OPEs but did 
not give the simple lagrangian description above. 

It is possible to construct "twist" fields |2§| for the field which play a role similar 
to the spin field of the Ising (or Majorana) field theory in the construction of the bulk 
quasihole wavefunctions J7[. These fields obey identical relations as those defined in the 
next subsection, so we postpone discussion until then. 



C. Relation of bulk and edge field theories 

Next we will explain briefly how the relation of the bulk and edge field theories can be 
used in order to define quasihole operators at the edge. 

The field theories as defined in §§IIIA, B appear very similar. The difference is that, in 
§IIIA the fields were not required to satisfy a reality condition, they and their adjoints both 
appeared in the action, and the Hilbert space was found to have positive norms (but no 
Lorentz invariance), while in §IIIB fields were real, Lorentz invariance was maintained and 
the self-overlaps of some states were negative. Here we will consider the positive definite 
theory of §IIIA, and exhibit a conformal structure in this system. This does not contradict 
the earlier statements because the stress-energy tensor involved is not self-adjoint (with 
respect to this inner product). 

The correlators of the (gradients of the) fields in this system, as already exhibited in 
( |3.12j ), are clearly conformally invariant. If we work in imaginary time r, and use the 
spacetime coordinate z = x + iv s r, then the fields obey the operator product expansion 
(ope) 
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dy a (z)d& a (0) ~ e aa ,/z 2 (3.15) 

up to the usual less singular terms (these have the same form as those for the real fields in the 
theory in §IIIB). Then, if we consider only correlators of d^ a , not o/cM/t, these correlations 
are conformally invariant. The stress-energy tensor that generates these transformations is 
T(z) = —e aa i : d^$! a d^ a ' : /2, which can be verified to obey the ope's of a stress-energy 
tensor using only the ope (|3.15|) . This operator is not self-adjoint, so its Fourier components 
L n in general do not obey L\ = L_ n . However, the Hamiltonian L derived from T coincides 
with that found above (|3.1| ). Naturally, if we consider instead dW a , there will be another 
non-hermitian stress tensor generating conformal transformations of those correlators. The 
vacuum is annihilated by the modes L n , n > — 1 of either of these stress tensors, as required 
in a conformal field theory. 

We may now consider twist operators in this theory, which quite generally are operators 
that twist the boundary conditions on the fields, in the manner already described in §11. We 
will require these to be also Virasoro primary conformal fields for the stress tensor T above, 
as they are in the non-positive definite theory. They are introduced by the operator product 
expansion 

d* a {z)S{w) ~ —L=g_ a (w) (3.16) 

y/Z — W 

where 5 is a twist field and S„ is an excited twist field of spin a. In the presence of the S 
fields the correlator for <9\I/j and which, in the untwisted sector, is given by the leading 
term in the ope ( |3.15| ) 

(0^(3)0^))= 1 (3.17) 
[z — w) z 



becomes 



(5(00)^(2)^(10)5(0)) = ' ( ^_ + ^' , (3.18) 



and by a standard calculation 03 we come to the equation 



(5(00)^)5(0)) = -l\ (3.19) 



which means that the conformal weight for the primary field 5 is —1/8, and that the corre- 
lator 

(5(*)5(0)) = z 1/4 (3.20) 

increases with separation. Negative scaling dimensions cannot appear in a conformal field 
theory described by unitary representations of the Virasoro algebra. Since, in theories 
obeying the BPZ axioms, operators correspond one-to-one with states, they indicate the 
existence of states with energy below that of the ground state. Since our Virasoro generators 
L n generate a nonunitary representation, this is not a problem here. Moreover, L coincides 
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with the physical Hamiltonian in both the untwisted and twisted sectors, at least up to an 
overall constant in the twisted case. (We note that in order to calculate this from AM of the 
zero energy states a term related to the contribution of the bulk to the angular momentum 
must be subtracted; see the next section.) We now propose that this constant is —2nv s /8L 
as predicted by the conformal considerations above, since this appears to contradict no 
principles. There is no real inconsistency in asserting that the energy of the supposedly 
"excited" twisted ground state lies below that of the untwisted "true" one. We could take 
the twisted sector (where ^ a obeys antiperiodic boundary conditions) as the ground state 
sector. We do not do so because (i) the theory for the disk (i.e. the chiral theory) clearly 
identifies the periodic sector as the ground state, which has no quasiholes in the bulk; (ii) 
the ground state in the antiperiodic sector is not invariant under SL(2, C) generated by L , 
L±i as required in a conformal theory, whereas that in the periodic sector is. 

For the edge states of the HR quantum Hall state (as in all the paired theories considered 
in this paper), the twisted sector of the fermions occurs only when the charge added at the 
edge is l/2q (modulo l/q). Making use of conformal arguments for the states with added 
charge, we expect the contribution to the energy from the charge sector be (2ttv c /L)Q 2 /2q, 
where Q/q is the charge added to the edge, as will be discussed in the next section. Thus, 
depending on the ratio v c /v s , the net energy of these states in the HR case will be positive 
in most cases, except for the lowest added charge sectors. In particular, for v = 1/2, the 
orginal case of interest for HR, there will be a sector with ground state energy — 2nv/16L 
if v s = v c . These sectors correspond to operators of the form Se lip ^ 2 ^, which are also the 
operators used in the bulk conformal field theory to generate quasihole states that are single 
valued with respect to the electron operators 0. 



IV. EDGE STATES AND THEIR FIELD THEORIES ON A CYLINDER 



A. General results and the Laughlin states 



In this section we consider zero energy states on a cylinder. For the Laughlin states, 
considered in the present subsection, the structure of the edge states is well known (see 
especially @), but will be reviewed here to ensure that the ideas are clear, and so as to 
introduce the partition function for two oppositely-moving edges. 

On a right cylinder of circumference L, we work in the Landau gauge. In terms of 
a complex coordinate z, the single particle wavefunctions in the lowest Landau level are 
e 2mnz/L e -^y ^ wnere n j s an integer and y = Imz; it has been assumed that the boundary 
condition is that wavefunctions are periodic under z — > z + L. A more general boundary 
condition is that the wavefunction changes by e 1 ^ under such a transformation; in that case 
the wavefunctions become e l ( 2 ' Kn +'t > ) z / L e -2y . Returning to the case e 1 ^ = from here on, we 
can write many-particle wavefunctions in terms of Zj = e 2mz ^ L , for example the Laughlin 



state 28 



i<j i 
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As any Z{ approaches any Zj, this function clearly retains the properties of the Laughlin 
state in the plane, namely it vanishes as the qth power. 

The Landau gauge has explicit symmetry under translation around the cylinder, and the 
corresponding conserved quantum number, M, is the sum of the powers Mj's of the Zj's 
in the many-particle wavefunctions. M is the angular momentum; alternatively we could 
use the linear momentum equal to 2nM/L. For e 1 ^ — 1, M is integral. In the Laughlin 
state above, M = qN(N — l)/2, and the Mj's of individual particles are in the range 0, 
. . .q(N — 1). The single particle wavefunction with Mj = n is peaked at y = ~2irn/L, so 
the Laughlin state occupies a corresponding range in the y direction. Due to translational 
symmetry in the y direction, there are, however, also an infinite number of other Laughlin 
states obtained by shifts, which are produced by acting repeatedly on the wavefunction with 
Hi Zi (or its inverse). This operator shifts the Mj by 1, so since the filling factor is 1/q, it 
corresponds to shifting a charge 1/q from one edge to the other. Fractional shifts change 
the boundary condition on the wavefunction, so are not allowed in the Hilbert space at fixed 

The infinite set of Laughlin states ("ground states") are the most dense or compact 
zero energy states, in the sense that the Mj's lie in a range of minimum possible width 
(note that the pseudopotential Hamiltonian in arbitrary geometry is defined in terms of the 
order of vanishing of the wavefunctions as two particles approach each other, which in the 
plane is equivalent to relative angular momentum). The total angular momentum of the 
Laughlin states is M = ^qN(N — 1) — pN, for the state where we have applied Y\Z~ P , p 
integral, to the Laughlin state (|4.1| ). The range of Mj values found in this state is then 
— p < Mi < q(N — 1) — p. (Note that for q(N — 1) even, we could choose p such that 
M = and obtain a state symmetrical about Mj = 0.) All other zero energy states have a 
broader range of Mj's, and are obtained as edge excitations of the two edges. In the N — > oo 
limit, the different Laughlin states become infinitely far apart in angular momentum, and 
the assignment of edge excitations as belonging to a particular ground state (from which 
their angular momentum differs by a finite amount) becomes unambiguous. There are then 
two sets of elementary edge excitations, out of which these excited states can be built, and 
as for the disk, these are linearly independent in the iV — > oo limit. The elementary bosons 
that create them are the operators s n = Z™ and s n = Z~ n {n > 0). We refer to the 
excitations created by the action of the s n 's, whose contribution to AM to t, the change in 
M relative to the corresponding Laughlin state, is positive, as right-moving. The other 
operators, s n 's, create excitations at the other edge, have AM to t < 0, and are viewed as 
left-moving. (Here right and left refer to the two directions parallel to the edge.) Thus, if 
we can split AM tot into AM — AM, where the two terms are the contributions at the two 
edges, then we would like to view AM + AM as the "pseudoenergy" (within a scale factor), 
where both terms are defined to be non-negative for the case in the present subsection. The 
direction of motion then follows from the group velocity, the change in pseudoenergy with 
momentum of an elementary excitation. 

A more realistic, and well-defined, method would be to introduce a Hamiltonian that 
breaks the degeneracy of the zero-energy ground states and edge excitations; for example, a 
parabolic confinement potential SjMf, would suffice, but unfortunately the wavefunctions 
that can be easily written down are not eigenstates of such a form, except for q = 1. We 
assume that as such a term is turned on and the eigenstates evolve, they stay in one-one 
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correspondence with those found here. We expect on physical grounds that excitations 
that move electrons further from the minimum of this potential, as the edge excitations 
do, have higher energy, and so for small \AM\ the bosons have a dispersion relation E ~ 
| AM | for either edge, and form the modes of a (non-chiral) scalar Bose field. The effect of 
the Hamiltonian on the low-lying states of the system can then be determined through a 
renormalization group analysis, as has been done for the present case already |J. In fact, 
the analysis of the effective field theory of low-lying excitations, and its operator content, 
given here is the basis for such an analysis for the paired states. 

If the electron number in the Laughlin ground state is changed by 1, the width changes 
by q units. We can view this as a new operation that adds charge to a single edge, without 
disturbing the bulk ground state. The full set of possible ground state systems can then be 
parametrized by the positions of the two edges, which are almost independent. Each edge 
can be shifted by q units without affecting the other. Shifts by m = 1, . . . , q — 1 must be 
performed on both edges together. A shift by q units at both edges is equivalent to removing 
an electron from one edge and inserting it in the ground state at the other. If we extend the 
idea of a charge sector to include in a single sector all those that differ by integral charges, 
then there are only q sectors, of charge 0, 1/q, 1 — 1/q (modulo integers), where the 
charge is shifted by the stated amount from one edge to the other. For the pseudopotential 
Hamiltonian without the confining potential, these q sectors are on an equal footing and the 
choice of zero is arbitrary. Above these ground states, bosonic excitations can be created at 
either edge, and have the same character in all sectors. Thus, the Hilbert space of the edge 
excitations can be written in the form 

q-l _ 

V = 0V r/ ,®V r/ , (4.2) 

r=0 

Each Hilbert space V r / g (V r / 9 ), r = 0, 1, . . . q — 1, is the span of the full set of states in 
the extended charge sectors at the right- (left-) moving edge. The conformal field theory 
of the edge states that includes the operators of charge 1 in the chiral algebra is known as 
the "rational torus" (see, e.g., 0). Loosely, the chiral algebra is the algebra of operators 
that affect only a single edge. There is such an algebra for both edges, the two algebras 
are isomorphic, and operators in one commute or anticommute with those in the other. 
Each Hilbert space V r / 9 (V r / 9 ), r = 0, 1, ... q — 1, is an irreducible representation of the 
fully-extended right- (left-) moving chiral algebra. 

To extend our definition of the pseudoenergy, which was given for the edge excitations 
of a single Laughlin state at fixed (but large) particle number, to the full set of sectors 
just described, we introduce an arbitrary reference ground state with N = N Q particles and 
shift p = po, so that the angular momentum M = ^qN (No — 1) — p N . We calculate 
AMtot for Laughlin states where an integral amount of charge has been added to a single 
edge, by changing N and adjusting the shift p from p such that either the maximum 
or minimum occupied Mj is unchanged compared with the reference state in the case of 
charge added to the left- or right-moving edge, respectively. After subtracting a quantity 
related to the bulk of the system, a step analogous to measuring momentum from the Fermi 
wavevector in a Fermi gas, this gives a formula for AM or AM, valid in these special cases. 
This will then be used in all the sectors. The occupied states in general lie in the interval 
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— p < Mi < qN(N — 1)/2— p. To add charge at the right-moving edge, we let A" = Nq + AN, 
and p = Pq. Then we calculate 

AM tot = qAN 2 /2 + AN[q(N - 1/2) - p }. (4.3) 

Therefore if we define E = AM tot — [q(N — 1/2) — p ] AN for the excitation pseudoenergy 
at the right edge, we obtain E = qAN 2 /2, for the Laughlin states, and this is consistent 
with result for the charge fluctuation bosons, which do not change the charge at the edge. 
Similarly, for the left edge, where we must use the adjusted shift p = p + qAN so that the 
right edge is unmoved, we obtain 

AM tot = -qAN 2 /2 + AN[-q/2 - p ) . (4.4) 

Then for the excitation pseudoenergy E at this edge we must use E = —AM — [q/2+po]AN. 
For each edge, the coefficient of iV is the mean of the angular momenta of the highest 
(respectively, lowest) occupied single particle states in the reference state with N and with 
Nq + 1 electrons, so it resembles the Fermi wavevector in a Fermi sea, on two sides of the 
Fermi sphere. We now use these formulas also for the ground (and excited) states in any 
charge sector, having a combination of electrons added to Nq, shifts of charge from one 
edge to the other, and charge fluctuations; for such states we can always calculate, for each 
edge separately, the amount of charge effectively added, which may now be fractional, but 
is always a multiple of \/q. 

All information about the number and quantum numbers of the edge excitations in the 
thermodynamic limit (taken with L 2 /N fixed) can be conveniently summarized in a partition 
function analogous to that in §IIA. The partition function is a double series in a complex 
parameter x and its complex conjugate x, which contains information about right and left 
movers, respectively. In fact, if H is the Hamiltonian whose eigenvalues are 2ttv/L times 
the sum of an E and an E found above, where v is the speed of propagation of the edge 
excitations, and taking x = x = e~ 2lT/3v ^ L , then this partition function is the Gibbs grand 
canonical partition function of statistical mechanics, Tre _/3H . (This is consistent with the 
assignments in §11 and §111.) The E's of the charge sectors we have found agree with the 
conformal weights that are found in the conformal field theory . The result is expected to 
be the same for a more realistic Hamiltonian with a confining potential, as discussed above. 
The partition function is defined as 

Z(x,x) = Trx E x E . (4.5) 

The structure of V given above now allows us to express Z in terms of the trace over each 
space ~V r / q , and we define (using the Euler partition sum p(n) for the bosons from §IIA) 

oo oo 
£ /q {x)= £ (± l ) m J H +r )^J| (l _ s n r i (46) 

m=— oo n=l 

and the complex conjugate for V r / g . The traces, such as the xt/q s i over ^ e right- and 
left-moving spaces are known as characters, since they are essentially the characters, in the 
algebraic sense, of the irreducible representations V r / g of the chiral algebra. (The greater 
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generality afforded by the insertion of ±1, and by allowing r to be an arbitrary real number, 
will be useful later.) The partition function can then be written as 

Z(*,*) = X>^(*)| 2 . (4.7) 

r=0 

Similar structures to those found here for the Laughlin states on a cylinder will be found 
for the other states in the following subsections; unfortunately, for the paired states, while 
the method for calculating E and E produces a similar contribution for the different charge 
sectors, it is too crude to produce the analogous energies that we expect to originate in the 
fermion sectors, such as the —1/8 discussed in the preceding section, though we believe that 
they could in principle be obtained in a refined calculation. 

B. Pfaffian state 

As for the Laughlin states, the zero-energy states on the cylinder can be obtained from 
those in §IIB by replacing Zi by Z iy the gaussian factor by exp[— | Yli vV\-> an d recalling that 
the exponents of the Zi's run over all the integers, positive and negative. The formula for 
AM of the edge excitations in the untwisted sector containing fermions only still applies, 

F 

AM tot = ^(n fe + l/2) (4.8) 

k=l 

but now nk > describes right-moving fermions, nu < left-moving. For the latter we can 
define AM = -AM tot . 

By inspection of the resulting states, we deduce that, in the untwisted sector, if N even is 
fixed at Nq, then the total number of fermions excited is even, and the parity of the number 
at each edge (i.e., whether it is even or odd) must be the same. But if we increase iV by 1, 
we must create or destroy a fermion at one edge, as well as increase the charge by 1, which 
can be done at the same edge without affecting the other. So the chiral algebra includes 
the operator ipe l ^ v which does this p]. Then the parity of the fermion numbers at the two 
edges can be opposite provided the parity of the integral amount of charge added relative to 
the reference state is also opposite. Applying the operation, or its adjoint, once more to the 
same edge, we find states where the charge has changed by or 2 relative to the reference 
state, but the number of fermions has the same parity, still without affecting the other edge. 
Thus all these states lie in the same extended charge sector, and similar results hold in the 
other untwisted sectors where a fractional charge m/q has been shifted from one edge to the 
other. A total of 2q sectors results from these considerations. 

The twisted sector is obtained from the untwisted by inserting a factor 

(N-F)/2 

II ( Z <r(F+2l-l) + Za(F+2l)) 
1=1 

in the states with F unpaired fermions (in the notation of §IIB) to produce the effective 
shift by a half-unit that transfers charge l/2q from one edge to the other. The angular 
momentum of the excited fermions is then 
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F 

AMtot = $> fc (4-9) 

k=l 

where once again can run over the integers. This time n& > represents right-movers, 
rt& < left-movers, n& = is the zero mode, which cannot be assigned to the right or left 
moving sectors. On the other hand, the requirement that the total number of fermions must 
be even whenever the charge added to the reference state is even, and odd when it is odd, 
still applies. Thus we have the general projection rule that applies in all the sectors, twisted 
and untwisted, that the total number of fermions created in right-moving, left-moving or 
zero modes together must be equal to the parity of the total charge added to the system at 
the two edges together. In the twisted sectors, it can be satisfied by allowing either parity of 
right- and of left-moving fermions in all charge sectors, then choosing the occupation number 
(either or 1) of the zero mode to obey the condition. Consequently, the distinction between 
"even" and "odd" sectors, that existed among the untwisted sectors and was responsible for 
the factor 2 in the 2q sectors, no longer applies, and there are just q distinct sectors (or 
irreducible representations of the chiral algebra). The total number of sectors is therefore 
3q, which, in line with the general connection between bulk and edge states made in j7|, is 
the same as the number of zero energy ground states found in the toroidal geometry JT2 



Note that in this case the description of the chiral algebra as affecting only a single edge 
is not quite correct, because the right-moving operator if), and its left-moving analogue if), 
each contain a term that changes the occupation number of the zero mode. Nonetheless, in 
the Majorana field theory, these operators anticommute; similarly, the algebra of operators 
assigned to one edge does (anti-) commute with those assigned to the other, even in the 
twisted sector. 

The calculation of the partition function, which formalizes the above remarks, is conve- 
niently performed in terms of characters. The basic objects are characters for the states at 
one edge that differ in charge only by integers, and, in the untwisted sector, the parity of the 
charge difference from that in the lowest energy state is equal to the parity of the change in 
fermion number. Characters for the fermions alone will be useful; these are, for untwisted 
(antiperiodic) boundary conditions, 



i 

2 



n(i+x n+i/2 )+ n(i-^ n+i/2 

1=0 n=0 

oo oo 

X\{l + x n+1 ' 2 )-X[{l-x n+1 l 2 ) 



n=0 n=0 



(4.10) 



which are respectively for even and odd numbers of Majorana-Weyl (MW) fermions. The 
subscripts are the conformal weights of the corresponding primary fields, or the energies of 
the ground state in each sector. In the twisted sector (periodic boundary conditions), there 
is only a single nonvanishing character: 



oo 



n=l 



(The zero mode is omitted here, as it will be accounted for separately, as already explained.) 
The constant 1/16 in E for the twisted ground state here is the analogue of those in the 
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different charge sectors as derived in §IVA, where it was mentioned that we cannot at present 
derive this one directly from our zero-energy wavefunctions. 1/16 is the conformal weight 
of the corresponding operator a, the spin field, which twists the boundary condition on 
the Majorana fermion, like the twist field S discussed in §111. These three expressions are 
well-known as the Virasoro characters of the critical two-dimensional Ising model, as well 
as in other contexts. The characters of the chiral algebra relevant to the edge states of the 
Pfaffian are 

X^cvcn.untwistcdM = §Xo™ ' ( X ) (xt/q( X ) + Xr/q( X )) + \xf/^{ X ) (xt/ q ( X ) ~ Xr/q( X )) 
Xr^odd.untwistcdl^) — \x¥/-l{ X ) {xt/q( X ) + X r /q( X )) + ^Xo^ ( X ) {Xr/q( X ) ~ X r /q( X )) 
xfr+l/2)/< ? ,twisted( a; ) = Xl/16 ( X )X~(r+l/2)/q ( X ) ■ (4-12) 

It can easily be seen that these expressions are sums over states with the necessary constraints 
on the combinations of fermion and charge states included, apart from those that enter on 
combining right and left movers and zero modes. The partition function is, finally, 



rPf, 



9-1 

^ ' ^ |Xr/q,even. untwisted I 
r=0 



,Pf 



X ) + Xr/<?,odd,untwistcd( ;r ) + X(r+l/2)/g,twisted ( X ) ^ 



(4.13) 



C. HR state 



As for the Pfafnan state, the edge states of the HR state on a cylinder can be deduced 
almost immediately from the results for the disk. In the HR case, the untwisted sector is 
found now to contain zero modes, while the twisted sector does not. The zero mode, like 
the nonzero, right and left moving modes, can be occupied by a spin up or a spin down 
fermion, or both. From these states we can deduce the projection rule. It has the same form 
as for the Pfaffian, in all sectors. The projection rule requires even fermion number when 
no charge has been added, and that a fermion is created or destroyed whenever a unit of 
charge is added to a single edge, so that the chiral algebra includes an operator d^f^e 1 ^, 
very similarly to the Pfaffian. We can therefore write down the characters without further 
ado. The characters for the fermions are: 



Xo( x ) 
Xi( x ) 
X-i/s( x ) 
xt/s( x ) 



na+^v+iK 1 -*' 

_n=l n=l 

oo oo 

na+^-iK 1 -*") 2 

_n=l n=l 

oo oo 



r 



-X 



-1/ 



-1/8 



\{{l+x n+l l 2 ) 2 +Y[{l-x n+l l 2 f 

1=0 n=0 

oo oo 

f[(i + x n+1 / 2 ) 2 - x\{i-x n+1 ' 2 Y 



Ln=0 



n=0 



(4.14) 
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Note that the first two are the even and odd fermion number states in the untwisted sector, 
omitting the zero modes, while the last two are the same for the twisted sector, and we 
have used the negative conformal weight of the twist fields in writing the latter. In this case 
we have maintained the distinction between even and odd fermion numbers in the sector 
that contains the zero modes, so as to exhibit its fate explicitly. (If desired, the following 
approach can also be taken for the Pfafnan state, and the expressions already given for the 
partition function can be derived in this manner, verifying the argument, given in words in 
the preceding text, that there are only 3q sectors.) The characters for the chiral algebra of 
the HR state are: 



Ar/q,ev,untw 



» = iXo ( X ) (Xt/q( x ) + X r /q( X )) + \xt ( X ) (xt/q( X ) ~ X r/q ( X )) 
X™od ) untw( a = \Xl{x) (Xt/ q ( x ) + Xr/q( X )) + |Xo ( X ) {Xr/q( X ) ~ Xr/q( X )) 



A.(r+l/2)/g,ev,tw 



( x ) = \X-l/*{ x ) (xtr+l/2)/ q ( X ) + X {r +l 



/2)/5 



X 



+ \xt/&{ X ) (X( r+ l/2)/q( X ) - X( r +l/2)/q( X )) 



X(r+l/2)/q,od,tw( X ) 



1,.* 



2 X3/8i X ) [X( r +l/2)/q( X ) + X( r +l/2)/q ( X )) 

r+l/2)/q( X ) X( r+ l/ 2 )/q( X )) ■ 



+ \x\/^ x 



(4.15) 



We may now form the partition function, by combining the sectors subject to the rules 
already mentioned. In particular, in the untwisted sector where the zero mode occurs, we 
may combine right- and left- moving sectors of the same parity, in which case the zero mode 
may be either unoccupied or doubly occupied, or we may combine sectors of opposite parity 
if the zero mode is occupied once, which may be with either spin. Thus we find for the 
partition function: 

Z RR (x,x) = 
<?-i 



r=0 



q-1 

E 

r=0 



E 2 |Xr/<j,ev,untw( a ') + 2 Xr/g,od,untw i X ) + 2 (Xr/q,ev,untw{ X )X' 



^r/g,od,untw 



(X) 



+ Xr/< ? ,od,untw(- r )X r /^,ev,untw( ;r )) + X(r+l/2)/g,cv,tw( ;r ) 
Xr/q,ev,untw\ X ) ^ Xr/q,od,untw\ X ) + X(r+l/2)/q,ev,tw \ X ) 



+ 



X(r+l/2)/g,od,tw( X ) 



+ 



X(r-+l/2)/q,od,tw( ;r y 



(4.16) 



which shows that there are in fact 4q sectors. The untwisted characters have combined into 
simpler ones, similarly to the twisted Pfafnan state characters: 



(4.17) 



n=l 



which, however, appear twice in Z. The 4q sectors show that there are 4q primary fields of 
the chiral algebra in the system, which (in the notation of §111, except that the fields 
S are now the nonchiral fields that act on both right and left movers) are 1 (the identity 
operator), # T ^ i; Se^ + ^/ 2 ^, S aW e^ + ^ 2 ^ (S aW is the twist field excited in both left and 
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right sectors, and a, a =|, |), and these operators times additional factors e' 1 ^^^ which 
shift charge 1/q from one edge to the other. All of these fields are spin singlets, except 
S a ^ which transforms as spin \ <g> \ — © 1. The field ^\^\ is not strictly a primary field 
since it has weight zero, but this is not an important distinction here. All other fields are 
descendants of these, that is they can be obtained by acting with operators in the chiral 
algebras; as a particular example, the state created from the untwisted ground state with 
unoccupied zero modes by the zero mode of times a unit charge at the right moving edge 
is obtained from by acting with cM^e 1 ^^. Acting again, with a similar operator, leads 
us to the identity, so the representations are not irreducible, which is a peculiarity of this 
system. 



D. 331 state and the hierarchy and its generalizations 



The edge states and the partition function for the 331 state are, by now, easily obtained. 
There are zero modes in the twisted sector, as for the Pfaffian, but there are two types of 
fermions (particles and antiparticles) as in the HR state. Similar selection rules governing 
even and odd fermion numbers apply as in the other cases. Accordingly, the partition 
function can be written using the Weyl (or chiral Dirac) characters: 

f 
f 



xT y \ x ) 


1 

~~ 2 


oo 

n 

.71=0 


oo 

(i + x n+i / 2 f + na 

n=0 


- x n+1/2 




1 

~~ 2 


OO oo 
JJ(1 + X n+V2)2 

_n=0 n=0 


- x n+1/2 


Weyl / \ 
Xl/8 (X) 


_ 1 T l/8 


oo oo 

na+^+na 

_n=l n=l 


-x n f 


x7/f(x) 


= -x 1 / 8 
2 


oo oo 

na+^-iia 

_n=l n=l 


-x n f 



(4.18) 

71=1 

The characters entering the partition function are: 

\xT y \x) (Xr/q( X ) + Xr/q( X )) + \X^l\ X ) {Xr/q( X ) ~ Xr/q( X )) 

lx7;i\ x ) (xU x ) + x; /q ( x )) + ho cy \ x ) (xU x ) - 

\X\/S ( X ) (X(r+l/2)/q( X ) + X( r +l/2)/q( X )) 
+ ^Xg/8 ( X ) (X( r +l/2)/q( X ) ~ X( r +l/2) /q( X )) 
\X$/i ( X ) (X(r+l/2)/q( X ) + X( r +l/2)/q( X )) 
+ IX7;I\ X ) (xU m/ q( X ) - X^yqW) ■ (4.19) 

The partition function is 

Z 331 (x,x) = 

q-l 



v 331 (r) 

Ar/g,ev,untwv / 

v 331 / \ 

Ar/g,od,untwV / 

X(r+l/2)/q,ev,tw( X ) 



Xfr+l/2)/q,od,tw( X ) 



r=0 



( \xl 3 /q,cv,untw( X ) + Xr/ 9 ,od,untw ( X ) + 2 X(r+l/2)/g,ev,tw ( X ) + X(r+l/2)/q,od,tw ( X ) 



(4.20) 
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Again the twisted terms have combined to form a simpler character, 



X 3 (r+l/2)/q,tw( X ) — (Xl/16( X )) 2 X(t+l/2)/g( X )' (4-21) 

and there are two distinct sectors with this character. The equality of some characters of 
distinct sectors may also happen with the rational torus or Laughlin state characters, for 
which the characters obey X^/ q ( x ) — = * = ^^-r)/< ? ( x )- 

By bosonization ||, the Dirac (or Weyl) characters can be written, using the Jacobi 
triple product formula, in terms of characters for a chiral boson with q = 1 (summed over 
charge sectors): 

X^Xx)± X Y;j\x)= X ± 1 (x) (4.22) 



MW, 
/16\ 



x)) =4/2)/i(z) (4-23) 



and so the 331 partition functions can be written in the form of sums for two boson fields, 
which is described in detail below. As mentioned in §11, the bosonized description of the 
field theories is closely related to the description of the bulk wavefunctions as two-component 
generalizations of the Laughlin states. The E's for the various sectors can be obtained in 
that description by an argument similar to that given in §IVA for the Laughlin state. Since 
this includes a contribution from the fermions as well as from the charge degrees of freedom, 
this holds out some hope that a derivation in the "pairing" representation of wavefunctions, 
which might also be applicable to the Pfaffian and HR states, should exist. 

Here we will give without proof the general results for the hierarchy and its generaliza- 
tions, restricting ourselves, for simplicity, to the case where the matrix G below is positive 
definite, with the 331 state as a special case. Physically, this is the case where all modes at 
the same edge propagate in the same direction. The other case has been discussed recently 
in PJ30|. The bosonized field theory for the right-moving edge can be formulated in 



terms of chiral boson fields <p a , which have correlators (in imaginary time) 

(<p a (z)(pp(0)) = -5 af) hiz (4.24) 

where a, j3 run from 1 to n, and in the hierarchy n is here the number of levels. The 
a = 1 component could be taken to be the density fluctuation field we have used up to now. 
The others represent the internal, neutral degrees of freedom. In the composite fermion 



approach, n is the number of Landau levels for the fermions [13], and the density mode is 
usually taken to be the sum of the ip a . The chiral operators that are allowed to be used 
at this edge without affecting the other (which generate the chiral algebra) are of the form 
e iv aVa (z) ^ we uge sum mation convention). The vectors v, whose components in the basis 
labelled by a are v a , lie on a (Bravais) lattice A in n-dimensional Euclidean space, that is 
they take the form of integral linear combinations of n linearly independent vectors e a . Thus 
the chiral algebra can be generated by e ±ieaoi,pa ^ z \ a — 1, . . . , n. The scalar products of the 
e a are G a b = e a ■ e^, which defines the Gram matrix of A; G is positive definite here because 
we assumed Euclidian space. In the (generalized) hierarchy theory, G is a matrix of integers, 
and so A is an integral lattice. If v = v a e a then we have v • v' = v a G a bv' b . Since the basis 
labelled by a (which is not an integral basis) is orthonormal, we find that the conformal 
weight of the chiral operators is v 2 /2, and so is either integral or half-integral. 
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The possible shifts of charge or the other U(l) quantum numbers related to the compo- 
nents of if a from one edge to the other are described by similar operators that act on both 
edges simultaneously, as in the states considered earlier. The right-moving part of such an 
operator is of the form e %WaVoL ^ , where w is a vector in the dual (or reciprocal) lattice A* 
of A; the dual lattice is defined as the set of all vectors w such that w ■ v = an integer for 
all v G A. Clearly A is a sublattice of A*. The right-moving conformal weight (or ground 
state energy in the corresponding charge sector) of the operators e iWc,tpc, ^ z > is again w 2 /2, 
which is the same (modulo 1) as the statistical parameter 8/2ir of the quasiparticles in the 
bulk of the same state, which are also labelled by vectors in A* (the "excitation lattice" in 



the terminology of The lattice A (the "condensate lattice") labels the combinations of 
quasiparticles that make up the possible order parameters. The relation of these to changes 
in the charges at one edge was mentioned in |TJJ , and forms the basis for the results quoted 



here. If we view the lattices as additive groups, we find that A has index detG in A*, and 
so the quotient group A* /A is a finite group of det G elements. The extended charge sectors 
are labelled by the possible shifts modulo fields in the chiral algebra, that is, by the elements 
of A*/A. These can be described by a set of vectors wa, A — 1, ... , detG, in A*, one in 
each coset of A. The simplest case is the Laughlin states, where n — 1, G — (g), A = y/qZ, 
A* = Z/y/g, and A*/A = Z g , which is equivalent to the description in §IVA. Another simple 
case is the integral quantum Hall effect, in which for n Landau levels, G is (in a convenient 
basis) the n x n identity matrix, so A is the n-dimensional simple (hyper-) cubic lattice, 
det G — 1, and A* = A. Thus, in this case, all edge excitations are just electrons in the 
various Landau levels. 

For the 331 states, the matrix G in the basis natural for the ground state wavefunction 



in the form (|2 . 24 ) is 



G 331 = q+ ] q -]). (4.25) 

The resulting lattice A embodies the projection rules by being distinct from an orthogonal 
direct sum of one dimensional lattices. We note that there are det G = Aq extended charge 
sectors, as found earlier. Moreover, unlike the lattices for the hierarchy states, for the 
generalized hierarchy states, and the 331 state in particular, the sublattice (A*)- 1 of operators 
that are neutral is not the same as the sublattice A 1 - of A. The lowest- weight neutral operator 
that is in A* but not in A represents the Fermi field ip a of weight 1/2. Since it does not 
appear in A, it cannot be applied to one edge, but must be combined with a similar operator 
for the other edge, or with another Fermi field or a charged field at the same edge. This is a 
consequence of the projection rules that we saw using the fermionic form of wavef unctions; 
states differing only by one fermion at one edge do not both exist in the Hilbert space. In 
contrast, for the hierarchy, the projection rules place no restriction on the neutral operators 
that can be applied to one edge, since all vectors in (A*)- 1 are also in A -1 . 

We may now describe the partition function for the (generalized) hierarchy states, in 
the case where G is positive definite. The sectors are labelled by the det G cosets of A, 
and in each sector the chiral characters are sums over vectors in the coset, together with 
fluctuations in the <p a . That is, define 
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X a(x) = £ (wa+v)2/2 J] (1 - a; m )-". (4.26) 

veA m=l 

The partition function is simply 

detG 

Z\x,x)= £ \Xa(x)\ 2 . (4.27) 

For the 331 states, this agrees with that derived by bosonization from Z 331 . 

It is interesting that, in the theory of the edge states, the relation of composite boson and 
composite fermion approaches maps exactly onto the usual 1 + 1-dimensional bosonization 
(or its inverse, fermionization). 



E. Orbifolds, chiral superalgebras, and modular transformations 

For readers familiar with, or ready to learn about, CFT, we mention that the theories 
for the cylinder described in this section are examples of the construction known as "orb- 
ifolding". Definitions, results and examples of orbifolds can be found in Refs. P5|J3T|,|3^[| . In 



brief, the general algebraic definition of an orbifold involves starting with a rational CFT 
with a chiral algebra A on which some finite group Q acts as a symmetry. One then takes 
the subalgebra Aq that is invariant under Q as the new chiral algebra. The representations 
of A will be representations of Aq also, but will in general be reducible; each irreducible 
component transforms as an irreducible representation of Q. In addition, there will be new 
"twisted" representations of Aq that are not representations of A. The same operations 
are applied to the left-moving chiral algebra A and its representations. The (symmetric, 
diagonal) orbifold CFT then has a primary field for each representation of Aq, which at the 
same time are primary for the isomorphic left-moving algebra Ao. The rule for combining 
left- and right-moving representations is that all fields must be invariant under the simul- 
taneous action of Q on left and right movers, and untwisted (twisted) fields must combine 
with untwisted (twisted). 

The CFTs of all the paired states described in this section are examples of orbifolds with 
Q = 7i2- The algebra A for the Pfaffian is generated by the fields e ±l ^ v , ip, and contains the 
U(l) current algebra generated by dp and the Virasoro algebra for %p as subalgebras. The 
primary fields are e tr(p '^, r — 0, 1, .. . , q— 1 (we suppress the left-moving operators for now). 
Z 2 acts simultaneously on <p, ip by p — > p + ir^/q and tp — ► — tp. The algebra Ao is generated 
by ipe ±i ^ tp , and the primary fields are e %r<f> /v*, ipe trv ^, r = 0, 1, . . . , q — 1, which result 
from the splitting of the representations of A (the untwisted representations), together with 
the twisted representations ae 1 ^^^ 2 ^^, r — 0, 1, . . . , q— 1, which include the spin field a, 
which is the analogue for the Majorana fermion of the twist field discussed in §111. The states 
(or descendant fields) in these representations obey the "projection rules" found earlier, and 
the full description of the combination of left and right movers, and the resulting partition 
functions, can be done in agreement with the rules obtained from the wavefunctions in this 
section. Very similar descriptions work for the HR and 331 states. (For the 331 state, the 
orbifold we find is that where ip^ transforms by e m and ipi transforms by e~ l7T . These factors 
are both equal to —1, but the point is that the factors written describe the way the phase of 
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either field winds on taking it round the twist field; there is also an adjoint twist field around 
which they wind the reverse way. This behavior is required by the structure of the twisted 
states, which have definite pseudospin as well as charge quantum numbers. It is most easily 
understood in the bosonized representation ip = ipx, ip^ = e tip2 , in the same notation as in 
§IVD. Then the symmetry is (pi — > (pi + ir^/q, <p>2 —> y?2 + 7r. This orbifold leads back to 
the lattice described in §IVD. In particular, the spin (or twist) field for ip a is bosonized as 
e^ 2 / 2 , which must appear in combination with some other charged fields, as for the other 
orbifolds and the generalized hierarchy theories in §IIID.) In all cases, the rationale for the 
structure is that the electron (or other fundamental charged particle) is represented in the 
edge theory by a field like ipe^^^, which has fixed boundary conditions in all sectors, and 
all fields must be local with respect to it, just as in the bulk, all wavefunctions must be 
single valued functions of the electron coordinates j7|. 

Our description of the orbifolds glossed over one aspect of the systems discussed here. 
The usual definition of chiral algebras assumes that all fields in the chiral algebras (both 
A and Aq) have integral conformal weight. In our examples, tp and tp a that appear for the 
Pfaffian and 331 states have half-odd- integral weight, and for q odd, so does e Thus A 

is actually a chiral superalgebra in these cases 0, and so is Ao in some cases (and also for 
the algebra of the Laughin state for q odd, and the generalized hierarchy states whenever 
applicable to electrons). We emphasize that, for our purposes, a superalgebra is one where 
some fields have half-odd- integral conformal weight, rather than one where some of the 
relations are anticommutators instead of commutators. In fact, to describe electrons, which 
are fermions, rather than the quantum Hall effect of charged bosons, the chiral algebra is 
always a superalgebra except in the case of the HR states, due to the violation of the spin 
statistics theorem there as discussed earlier. 

The fact that the chiral algebra is sometimes a superalgebra has consequences for the 
modular transformation properties that we may expect for the partition functions calculated 
in this section. If x = e 2ntT , and Imr > (r should not be confused with earlier uses of the 
same symbol), then modular transformations act as 

ar + b , 

r ^ 4.28 

cr + d 



and the matrix 




(4.29) 



is a member of SL(2, Z), the group of 2 x 2 integer matrices of determinant 1. (The group 
of modular transformations themselves is SL(2, Z)/{±/}.) The modular group is generated 
by the elements T : r — > r + 1, represented by 



and S : r 




(4.30) 



(4.31) 
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When the chiral algebra is strictly an algebra (i.e. not a superalgebra) , then the parti- 
tion functions will be modular invariant, if we modify the definition to include the factor 
(xx)~ c * 2A ) where c is the central charge of the CFT (not the matrix element just above). 
Central charges are additive; the values of the central charge are c = 1 for the Laughlin 
state, 1 + 1/2 = 3/2 for the Pfaffian, 1 - 2 = -1 for the HR, 1 + 1 = 2 for the 331 (all 
independent of the value of q), and n for the (generalized) hierarchy states. Modular invari- 
ance occurs for q even in the case of the Laughlin and HR states, and q odd for the Pfaffian 
and 331 states, all of which except the HR state describe the fractional quantum Hall effect 
of charged bosons, not electrons. When some of the fields that generate the chiral algebra 
have half- integral conformal weight, they will obey an antiperiodic boundary condition in 
the space direction, and the (modified) partition function cannot be invariant under the 
full modular group; the only boundary conditions that are invariant under the whole group 
are periodic around any cycle on the torus. In these cases, which as we have seen apply 
to all states considered here that can describe electrons, with the sole exception of the HR 
states with q even, we expect that our expressions are invariant only under the subgroup 
of the modular group that leaves the antiperiodic boundary condition on the electron field 
invariant. This subgroup is generated by the elements S and T 2 and can be shown to be 
isomorphic to r (2)/{±J}, where T (2) is the subgroup of SL(2, Z) consisting of matrices 
where the matrix element c = (mod 2). 

We also mention here some isomorphisms of the chiral algebras of our systems to known 
algebras. (In this paragraph, c is the central charge, and N is not the number of particles.) 
For the Laughlin state (of bosons) at v = 1/2, the fields e^^ 2 ^, dip generate the SU(2) 
current (Kac- Moody) algebra of level 1. For the Laughlin state at v — 1/3, we have the 
N = 2 superconformal algebra at k = 1, generated by e ±iv ^, dip. For the Pfaffian state 
(of bosons) with q — 1, the operators e ±i v / ?v are the bosonized representation of a Dirac 
field, or of a pair of Majorana fields ip±i, which together with the Majorana field ip = ipo 
forms a triplet of Majoranas. This c = 3/2 theory contains an SU(2) current algebra of 
level 2, or equivalently an 0(3) algebra of level 1, in which the currents are the bilinears 
ipaipb, a, b = ±1, 0. This symmetry shows up, for example, in the degeneracies of the 
excited energy levels, as long as the velocities for ip and ip are equal. The 3q = 3 sectors, 
even untwisted, odd untwisted, and twisted, correspond to primary fields that transform 
respectively as spins 0, 1, 1/2, under both the left- and right-moving SU(2). Moreover, 
the product ipiipoip_i generates N = 1 superconformal symmetry ]3"3] ], though this operator 
does not survive the projection to A$. Finally, the algebra Aq for the is = 1/2 Pfaffian 
state is generated by ipe ±lvf2ip , which has weight 3/2, and the algebra can be recognized 
as superconformal iV = 2 at k = 2 |33]| . In this case, the unprojected algebra A contains 
SU(2) level 1 and an SU(2) triplet of supercurrents tpe^^ 2 ^, ipdip, which generate an N = 3 
superconformal algebra j|33|| . 



V. CONCLUSION 

To conclude, we have found complete descriptions of the wavefunctions, the Hilbert 
spaces and the field theories of the edge states of the paired systems considered. The explicit 
wavefunctions are very appealing and make the enumeration of excited states in terms of 
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elementary excitations straightforward. The combination of complete proofs of some results, 
and enumeration for low excited states in others, makes the correctness of those results not 
explicitly proven here almost certain. The results confirm the general prediction in J7| of a 
relation between bulk and edge properties. For the Pfaffian and HR states, this provides 
indirect evidence for the prediction of nonabelian statistics of the quasiparticles in the bulk 
of these states. The twist fields in the edge conformal field theories proposed here certainly 
have such properties when exchanged in spacetime at the edge ("monodromy"). For the 331 
state, as for all generalized hierarchy states, the monodromy, and the statistics of the bulk 
quasiparticles, is abelian. 

The explicit wavefunctions for the edge excitations, particularly for two edges on a cylin- 
der, are reminiscent of results for integrable one-dimensional systems, especially those of the 
Calogero-Sutherland (CS) type [E3] for which there are explicit, simple wavefunctions for the 



ground state and many excited energy eigenstates. Indeed, the similarity of the Laughlin 
state and the ground state of the CS model has often been remarked. In the limit L ^> N, 
the Laughlin state on the cylinder essentially becomes the Calogero-Sutherland ground state 
27J ; the coordinate transverse to the edges can be viewed as the canonical momentum. The 



edge excitations of the Laughlin state are in one-one correspondence with the excitations of 
the CS model, and the low-energy CFT of the latter is once again a Luttinger liquid. It is 
interesting to speculate that there might be some integrable one-dimensional Hamiltonians 
with long-range interactions, generalizing the CS model, for which the ground and excited 
states might be related in a similar way to the wavefunctions discussed in this paper. If 
so, then we expect the low-energy field theories of the one- dimensional models to be the Z2 
orbifolds discussed in §IV. 

Finally, we note that the approach used here can be applied to other states for which 
the ground state is the zero-energy eigenstate of a suitable local Hamiltonian, as here. An 
example is another paired state, the permanent state J7|, which is a spin singlet, and is the 
densest zero-energy eigenstate of a certain 3-body Hamiltonian []55fl . The resulting theory 
is a Z2 orbifold containing spin 1/2 bosons of conformal weight 1/2 at the edge. Such a 
system, like the HR state, violates the spin-statistics connection, so the edge field theory 
is not conformal. The corresponding non-positive conformal field theory, whose correlators 
reproduce the bulk wavefunctions, is the (3-^ ghost system J7|, so the relation of bulk and 
edge theories is maintained. 
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APPENDIX A: ZERO ENERGY STATES FOR THE THREE-BODY AND 

HOLLOW-CORE HAMILTONIANS 

In this Appendix we will justify directly the general form of the zero energy eigenstates 
of (|2.10| ), and by extension its analogues for q > 1, and show in particular that they lead 



to the forms for the edge states in ( |2.12p , ( |2.15| ). We then briefly address similar questions 
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for the hollow-core Hamiltonian for which the HR state is the unique ground state, and 
corresponding issues on the cylinder. 

The Hamiltonian ( E7T0 ), taken with Bose statistics for the particles so that the Pfaffian 
state with q = 1 is a possible ground state, implies that the wavefunction of a zero energy 
state vanishes whenever any three (or more) particles coincide. This implies that zero energy 
states can be written in the form of the Vandermonde determinant lli<j( 2 i ~~ z j)i times an 
antisymmetric function that, as a function of any two coordinates Zi, Zj, may have a simple 
pole at Zi = Zj, times the usual gaussian factors. Such a state will have zero energy provided 
the antisymmetric function involved does not have a triple pole of the form 

[(Zi - Zj)(Zj - z k ){z k - z^ 1 

as i, j, k approach one another, for any i, j, k. A form like 

[(Zi - Zj){Zi - Zk)}- 1 

cannot appear either, because it is symmetric in j and k, while a form 

(zj - z k )(zi - z j y l (z i - Zk)' 1 

could, but this can be rewritten as a difference of simple poles 

(Zi — Zj) 1 — [zi — z k ) 1 . 

So all possible functions can be written as linear combinations of the forms already given, 
where the singularities involve disjoint pairs of particles. 

Without loss of generality, the general state can be taken to be a linear combination of 
states written by antisymmetrizing a function obtained by dividing the particles into pairs, 
writing an odd factor for each pair and symmetrizing over exchange of pairs. These condi- 
tions are of course sufficient but not necessary for the final antisymmetrization over all parti- 
cles to be nonvanishing. That is, neglecting the omnipresent factor IK^i - z j) ex P(~ \ \ z i\ 2 )-, 
we must have 

v Y,raS N/2 rifeLl fk( z <r(2r(k)-l), z a(2r(k))) 

L s s na - ( — \ — 7 r ( A1 ) 

a&S N 1^(1) - Z *(2)J ■ ■ • \MN-1) - z a(N)) 

where the f k are symmetric polynomials in two variables. For iV odd we can write a similar 
form with k — 1, . . . , (N—l) /2, and include for the unpaired particle an arbitrary polynomial 
factor f (z a{N) ). 

A convenient way to describe the symmetric functions in two variables z±, z 2 , is the 
following. We know that symmetric functions can be written as sums of products of the 
sums of powers s$ = z™ + z%- We will separate the symmetric functions that vanish at 
z\ = z 2 by writing the disjoint sets of functions 

A m = {(*i - ^) 2m 4 2) = n = 0, 1, 2, . . .} (A2) 

for m = 0, 1, 2, .... We claim that the full set of symmetric functions in two variables is 
spanned by linear combinations of the polynomials in the set 
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(J An (A3) 



m=0 



(there is no need to take products of these functions). This can be shown by induction from 
the fact that products of sums of powers span the symmetric polynomials, together with the 
identities 



sgag = 2s% +n2 - {zT - z^){zT - zT) (A4) 
( Zl - z 2 ){zr l + z^ 2 z 2 + -.. + zT 1 ) (A5) 



which (by induction on the order) express a product of elements of A as a linear combination 

of elements of A , A\, 

Now each fk in (|A~T| ) can be chosen to be an element of Um=o An- If fk is an element 



sffl of A , we will try to pull outside the sum on permutations a the corresponding sum 
of powers in all A" coordinates, s n (see eq. ( |2.6| )); this will leave behind terms with fewer 
fk G Aq. Repeating this procedure, eventually all fkS remaining inside the sum will be in 
Um=i A m and we will have finished, fk that are in A m {m > 1) contain {z a {2 T {h)-i) — z a (2 r (k))) 2 
which cancels a factor in the denominator, so these particles are unpaired in this term and 
the wavefunction will be a linear combination of the forms (|2.12|) . 

For the basic (untwisted) sector of edge states, we can consider A" large and most fk = 1, 
though this is not necessary and the results below are valid for all wavefunctions of the 
stated form. Then we observe that if /i = s%> G A , then 

N/2 N/2 N/2 N/2 

e ru«*n e ru- e ea n /* cab) 

tGS n / 2 fc=l TeS N /2 k=2 reS N / 2 k=2 k'=2,k'^k 



where for k = 2, 3, ... , 



f - ' l{fk = l (A7) 
fk ~ 1 /,4 2) if /* ^ 1, ( 7) 



and s n can be taken outside the sum on permutations a. The functions fksffi can then be 
reduced using the identity ( |A4|) , and all the terms in the many-particle state are now of the 
form of symmetric polynomials in A^ variables times antisymmetric functions with fewer fk 
that are members of Aq and ^ 1. Eventually, all fk are either 1 or are G A m (m > 1), and 
these states are linear combinations of the states in the text. Similar methods work for A^ 
odd. Thus we have shown that all zero energy states are linear combinations of symmetric 
polynomials times the form in ( 2.12j) . 



To obtain the twisted sector we can replace fk by fkSi in the above proof, leaving the 
s\ factors intact inside the sum on a and r at each step. Of course, in a finite system, our 
proof shows that these states can be expressed as combinations of the others, but to study 
the Hilbert spaces of edge states, we take A^ — >• oo before the number of fk ^ 1 becomes 
large, and thus we obtain two different sectors in this limit. Similar arguments apply if it is 
desired to include any other factor in every fk in the state. 



The hollow-core Hamiltonian [10] requires that zero energy states have no pairs of par- 



ticles with relative angular momentum q — 1. (Another way to say this, which is useful in 
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other geometries, is in terms of the order of vanishing of the functions.) We recall that the 
relative angular momentum of a pair, say 1, 2, is defined by expressing the wavefunction in 
the form (neglecting the gaussian factor, and the spin labels if any) 

oo 

z 2 , z 3 ,...,z N ) = (*i - z 2 ) m (z 1 + z 2 ) n ^ mn (z 3 , ...,z N ) (A8) 

771,71=0 

in which each term in the sum is an eigenstate of relative angular momentum of 1 and 2 
of eigenvalue m. The densest zero energy state of the hollow-core Hamiltonian occurs at 
filling factor v = 1/q. The largest q for which the pairing in which we are interested can 
occur is q = 2. For q > 2, zero energy states can be obtained from those for q = 2 by 
multiplying by the Vandermonde determinant. For q = 2, the wavefunctions are required 
to be totally antisymmetric when the spin states are included (see §IIC). For fixed spins of 
the N particles, the wavefunctions can be written as Yii<j( z i ~ z j) 2 times a meromorphic 
function, and the meromorphic function must be antisymmetric among particles of the same 
spin. For zero energy states this function must have, as any two particles come to the same 
point, either a double pole, with zero residue, or be analytic. Because of antisymmetry, 
double poles can appear only for opposite spin particles. All antisymmetric functions can 
be obtained by antisymmetrization of functions of indefinite symmetry, though we may as 
well omit functions that would vanish on antisymmetrization. If a double pole is present 
for a pair i, j, then it cannot be present for any other pair of the form i, k or j, k. This is 
because the Vandermonde squared contains the factors 



2 



2 



(zi - Zj) 2 (zi - z k ) 2 (zj - z k f = (zi - zj) 2 I y\(zi + Zj) - z k \ -\{zi- Zj) 2 j (A9) 

the expansion of which contributes only even powers to the relative angular momentum of % 
and j. In view of the double pole in Zi—Zj, the whole wavefunction is a zero energy eigenstate, 
provided there is not a single or double pole in Zj — z k or Zj — z k , for any k. Therefore, all 
pairing factors (zj — Zj)~ 2 must contain distinct pairs. The unantisymmetrized function can 
thus be written as a product of pair factors for as many opposite spin pairs as possible, 
times functions fk( z J, z j) of the paired coordinates that can be taken either symmetric or 
antisymmetric, and must be either nonvanishing at zj = zj, or vanish at least as fast as 
(zj — z^) 2 , so as not to spoil the zero energy property. The zero energy wavefunctions thus 
have, without loss of generality, a form similar to (IA1I). We now try to pull any one of the f k 

T I 

that does not vanish at z\ = Zj outside the sum on permutations by the same procedure as 
for the Pfaffian, using ( |A6| ) . Use of (A3), ( |A5| ) then shows that the resulting functions still 



obey the zero energy property. The procedure can then be repeated until linear combinations 
of the form (|2.20|) are reached. We conclude that the wavefunctions in the form ([2.20 ) span 



all the zero energy states, in the untwisted sector. Similar arguments apply to the twisted 
sector, to 7^ iVj, and to combinations of these. 

Finally, we comment on zero energy states on the cylinder. On replacing z^ by Zi (see 
§IV) we see that in (|A1|) f k must still be symmetric but may contain negative powers of 
Zi. (The pairing factors (Zi — Zj)' 1 can be left unchanged without loss of generality.) We 
extend the definition of the sets of symmetric functions (|A2|) by allowing the exponents n 
in the symmetric polynomials sffi to be negative as well as positive or zero, while m is still 
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non-negative; we claim that these span all symmetric holomorphic functions in two variables 
on the cylinder. (|A4j) , (|A5[) apply unchanged to all integral values of n%, n 2 , although (|A5|) 
becomes an infinite series. The proof then works as before, by pulling sums of (positive or 
negative) powers outside the sum on permutations. The HR case works similarly. 



APPENDIX B: LINEAR INDEPENDENCE FOR SMALL AM IN THE PFAFFIAN 

AND HR CASES 



In this appendix we construct and verify the linear independence of the states in the 
untwisted, even N sector for all AM < 8 for the Pfaffian and AM < 6 for the HR state. We 
use a different basis from that derived in Appendix A. We first return to the two-quasihole 
states for the Pfaffian. Due to the symmetry of exchanging w± and w 2 , they may be expanded 
in the form: 

N/2 m 

<if(z 1 ,...,Z N ;W 1 ,W 2 )= Y l T / ^rnn(Zl,...,Z N )(w^- n + W^- n ) (Bl) 

m=0 n=0 

where all the \l/ mn are linearly independent. This may be interpreted as saying that the 
quasiholes behave as two bosons, which may each occupy any one of N/2+1 states. This does 
not, however, mean that the quasiholes are bosons in general, which would contradict the 



assertion that they obey nonabelian statistics [35|. To obtain the expansion, we first expand 



the numerator in ( |2.11[ ) inside the Pfaffian, i.e. for a fixed choice of pairs, described by the 
permutation a (each pairing is obtained from 2 N / 2 (N/2)\ different a's). For each pair a(2k — 
1), a (2k) the factor [(z a{2k -i) - Wi)(z a ( 2k ) ~ w 2 ) + (wi <-> w 2 )] will contribute z a(2k _ i)Z a{2k) , 
(z a (2k~\) + ^o-(2fc)) or a constant to an expansion coefficient \I/ mn . This observation suggests 
the use of an alternative basis for the space of edge states spanned by the \l/ mn , defined by 



$AM,s(^l, ...,Z N ) = 

1 x - sgna 

2^(N/2)\ ^ (z a{1) - z a{2) ) ■ ■ ■ (z CT (jv-i) - z a{N) ) U 



• • • (z CT (iV-l)^(iV)) mjV/2 (^(l) + Z*(2)) ni ■ ■ ■ (Za(N-l) + ^(iV))™^ 2 }} 

xn(^-^exp[-V|^| 2 ] (B2) 

i<j 

where AM = M — Mq is again the difference between the total angular momentum M of 
the edge state and the angular momentum of the ground state Mq. The expression in the 
double curly brackets is defined as the sum over permutations of N/2 pairs: 

{{(^(l)^(2)) mi • • • (Za(N-l)Z a( N)) mN "(z a(l) + Z a(2) )^ ■ ■ ■ (z CT (AT-l) + Z a(N) ) n W } } = 

AT 1 J2 (z a(1) Z a{2) ) m ^) ■ ■ ■ (z a{N _ 1)Z(7{N) ) m ^ (Z a{1) + Z a{2) ) n ^) ■ ■ ■ (Z^.!) + Z a{N) ) n ^ 



tSSjv/2 



(B3) 



which makes this expression invariant under permuations of the pairs, and under permu- 
tations n a , m a 1— > n T '( a ), "V(a)- N is the number of permutations in Sn/2 that leave 
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the sequence of pairs n a , m a , a = 1, . . . , N/2 invariant. In the states $am,s the num- 
bers n a , m a are defined to be or 1, such that n a + m a < 1, and J2a=i m a = AM — s, 
YlH=i n a = 2s — AM. With these restrctions there is clearly just one distinct polynomial 
of the form (|B3|) for each s, which will be denoted Pam,s, and we see that s < AM < 2s, 
AM - s < N/2, 2s - AM < N/2, and s < N/2. Comparing $am,s and # mn , we see that 
AM = N — m. 

From (B2) it is easy to calculate how many edge states of fixed AM the expansion of 
two quasiholes gives, as N — > oo. There are 1 + AM/2 linearly independent states for AM 
even and (AM + l)/2 for AM odd. But for fixed AM 

AM 

p AA/, s = e AM (B4) 

s>AM/2 

where cam is an elementary symmetric polynomial, independent of the permutation a, which 
can therefore be brought outside the sum on permutations as a multiplicative factor. This 
arises because bringing the two quasiholes to the same position, W\ = w 2 , produces a single 
Laughlin quasihole. The remaining edge states, which span spaces of dimensions AM/2 
for AM even and (AM — l)/2 for AM odd, require nontrivial factors inside the sum over 
permutations and represent degrees of freedom that are not simply density fluctuations at 
the edge. At each AM, the number of such states coincides with the number of states with 
two fermions added to the ground state and the same AM in the Majorana field theory. 

Further expansions of the states with more than two quasiholes should generate, besides 
further symmetric polynomial factors, all even-fermion-number excitations. Wen I?| has 



demonstrated numerically for (|2.10|) for up to N = 10 particles that the number of low AM 



zero energy states coincides with that in the Majorana field theory. We will go a little further 
analytically, for arbitrary N. States with 2n quasiholes, n > 1, at positions {wi, . . . , w 2n } 
will, in place of the factor 

N/2 

II [( 2 ^(2fc-i) - Wi){z a[ 2k) - w 2 ) + (wi <-> w 2 )} (B5) 

k=l 



inside the sum on permutations a in ( |2.11| ), have a product of such factors, each involving a 



distinct pair of u>'s. Thus the degree of the wavefunction will be N<j> = q(N — 1) — 1+n. There 
are 2 n n\ distinct ways to associate the w's in pairs, but only 2 n ~ 1 of the resulting electron 



wavefunctions are linearly independent (as functions of the z^s for fixed w's) for n > 1 |p5 
Provided that they are degenerate in energy, which is true by inspection for the appropriate 
3-body Hamiltonian, the fact that this number is > 1 is the basis for nonabelian statistics. 
When the quasiholes are exchanged adiabatically, the usual Berry phase is replaced by a 
matrix acting in this space of degenerate quasihole states; however, this has not yet been 
explicitly demonstrated in this or any other example (see |^|). Here we are interested to 
see what edge excitations we can obtain by expanding these states. Evidently expanding 
in powers of Wi, w 2 n will generate the general polynomials (|B3| ) inside the sum on 
permutations, without restrictions on the m a 's and n a 's. We will then have to take into 
account the linear relations among some of these states, corresponding to those discussed in 
for n = 2. 
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The linear relations among some of the states are obtained from the following general 
identity. For any set of complex numbers a^, % — 1, . . ., P, P > 2 even, 



Pf(a i -o i ) = 0. (B6) 

This follows because the Pfaffian is the square root of a determinant in which any three rows 
or columns obey a linear relation. All cases with P > 4 can be viewed as applications of the 
identity for P = 4. One consequence of the identity is that when we insert the expression 

{{(Ml) - ^(2)) 2 (^(3) - 2a(4)) 2 }} (B7) 

(defined in analogy with (|B3|) ), or similar expressions, into the Pfaffian the resulting ex- 
pression vanishes. ( |B~7| ) can be expanded in the basis ( |B3| ) and this gives a linear relation 
among the states obtained. Therefore for AM = 4 we "lose" one state. For AM = 5 the 
expressions are 

{{{z a{1) - z ai2 )) 2 (z a{z) - z aW ) 2 (z a{5) + z ai6) )}} (B8) 

{{(z ff( i) - z a{2) ) 2 (z a{3) - ^ (4 )) 2 }}ei. (B9) 

Notice that in the first expression the distinct pair er(5), er(6) is introduced, which does 
not affect the vanishing which is due to the summation over permutations of the other four 
particles. In the second we have simply multiplied the insertion (|B7|) by an elementary 
symmetric polynomial, which is linearly independent of the other expression. From here 
on we will omit the expressions which are products of the ones valid for lower AM and 
symmetric polynomials. Then for AM = 6 the linearly independent expressions producing 
linear relations among the states are: 

{{(z a{1) - z a{2) ) 2 ( Z(T{3) - z ai4) ) 2 (z a(5) - z a{6) ) 2 }}, (BIO) 

{{(Z a{1) - Z a(2) ) 2 (z a(3) - Z a{4) ) 2 (z a{5) + ^(6)) 2 }}, (Bll) 



~~ z <r(2)) 2 (M3) - ^(4)) 2 (2 ff (5) + Za(e))( Z <r(7) + MB))}} ( B12 ) 



and 



UOa(l) - ^(2)) 2 (^(1) + M2))(^(3) - ^7(4) ) 2 (^(3) + ^(4))}}- (B13) 

The last expression uses ( |B6| ) with aj = zf. For AM = 7: 



- ^(2)) 2 (^(3) - ^(4)) 2 (^a(5) ~ ^a(6)) 2 (^a(7) + 2<t(8))}}, ( B14 ) 
{{{ z tr(l) - 2^(2)) 2 (^7(3) - 2<r(4)) 2 (z«r(5) + ^(6)) 3 }}, (B15) 
{{(Ml) ~ Z «(2)) 2 ( Z v(3) ~ ^(4)) 2 (^(5) + Z*(6)) 2 (Za(7) + ^(8))}}> ( B1 6) 
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~ ^(2)) 2 (2<7(3) - ^(4)) 2 (^ CT (5) + ^(6))(^(7) + ^(8))(^a(9) + ^a(lO))}} (B17) 



and 

{{(^(1) - ^t(2)) 2 0<t(1) + ^(2))(^(3) - ^a(4)) 2 (^ CT (3) + ^t(4))(^t(5) + ^(6))}}- (B18) 

All these expressions are linearly independent of each other. 

Since symmetric polynomials can always be multiplied into zero energy states to ob- 
tain another zero energy state, it is convenient to decompose all states into a product of a 
symmetric polynomial and another part that is linearly independent of symmetric polyno- 
mials. The latter represents excitations that are not density fluctuations at the edge. The 
full Hilbert space of edge excitations thus can be written as a tensor product of a bosonic 
Fock space of density excitations, as described earlier, and another space of independent 
excitations. Since the AMs of the excitations add, the dimension of the full space at any 
AM can be obtained by convoluting those of the two factor spaces. It is easy to calculate 
the dimensions obtained for the latter space by building its states up from products of the 
Pam,s (rendered linearly independent of e^f) an d then subtracting the number of linear 
relations just obtained. For AM < 7, we find that the linear relations eliminate all the 
states obtained from more than two quasiholes. Thus we find that for AM < 7 the edge 
excitations of the Pfaffian state exactly match those in the chiral boson times Majorana 
fermion system, in the fermion number zero or two sectors. 

In principle, it is possible to find the number of the edge states at arbitrarily high AM, 
by deriving these expressions in a systematic way. First we list all polynomials of degree 
AM of the form 

{{0<x(i) - z<7(2)) 2 0a(3) - M4)) 2 ■ • ■}} (B19) 

where dots denote additional squared differences or sums that multiply the first two terms. 
Then we take the space of expressions that vanish when inserted in the Pfaffian at lower 
AM, multiplied with all possible products of symmetric polynomials that make the degree 
of the expression AM. We expand these in the terms of the form ( [B19[ ). Some terms in 



the expansion are obviously zero when inserted in the Pfaffian; then the rest must give zero 
too. Fortunately for low momenta (AM < 7) each wave function of the form (P19|) is zero 
and that leaves us to prove only that the rest of the wavefunctions are non-zero and linearly 
independent. This can be done by taking pairwise limits z\ — ► Z2, etc, of particle coordinates 
in the Pfaffian alone (i.e. without Laughlin-Jastrow factor), whenever these are singular, and 
examining the linear independence of the resulting functions of the remaining variables, that 
are the residues of these poles. 

We will just state that the number of the edge states we found at AM = 8 implies that 
the four fermion state (|, l|, 2~, 3^) is present in the spectrum. So the numbers found are 
those given in the table in §IIB. 

We now turn to the HR state. States with any even number, 2n, of quasiholes are zero 
energy states and can serve as generating functions for edge states. As in the Pfaffian case 
they suggest an overcomplete basis of the states in which polynomials of the ( p3|) type 
are inserted. For more than two quasiholes, linear dependences arise when more than one 
factor like (zj — zj) 2 cancels a similar factor in the denominator; this time, the factors in 
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the denominator are themselves squared, and the sum over permutations of the relevant 
particles gives a determinant, not a Pfaffian, so the identity that replaces (B6) is simply 
that 



det 



/I 1 
1 1 



1\ 



1/ 



(B20) 



The linear dependencies that we mentioned in the Pfaffian case with (i f a(i) J,) pairs instead 
of (a(i)a(i + 1)) are then valid here, except the last one in the AM = 6 and AM = 7 cases. 
Namely, if we insert any of these in the HR state fl2.18|) the sum over permutations will 
produce zero. Therefore, the numbers of nontrivial edge states, that is, the states without 
symmetric polynomials we found so far in the Haldane-Rezayi case are 



AM 


1 


2 


3 


4 


5 


6 


7 


dim 





1 


1 


2 


2 


4 


4 



All these states are singlets because they consist of singlet pairs. 

The complete spectrum of edge states of the HR system should contain also nonzero 
spin states. At present we have no quasihole-like generating functions for these, but we can 
still obtain the edge states by writing down suitable functions directly. Each such state is 
formed when some of spin-singlet pairs in the ground state are excited into triplet states 
of zero energy. The broadest class of edge states made in this way, with S z = 0, has this 
polynomial in the numerator of the term with a fixed arrangement of coordinates in pairs: 

{{(4 - 4m) ni • ■ ■ (4/2 - 4(W njv/2 (4 + 4m) mi ■ • • (4/2 + 4 W2 )) mjv/2 }} (B21) 



where no rii = 1, and at least one of n'^s is an odd number. 

For AM = 3 we have only {{(zj — ^(i)) 3 }} which corresponds to S — 1,S Z = 0, i.e., 
belongs to the permutation group representation (2 N ^ 2 ~ 1 1 2 ). Similarly the two remaining 
states of the triplet can be constructed. For example the S z — — 1 state is: 

tf s ,=-i(4---,4/2 + i) = 
E s g ncr (4(i) 



reSjv/2+i 



Kj 



M2) \z\- 
Zi - Zj) q exp[-~J2\ z i 



M3)> 



- 7 l ) 2 
N/2-1 <t(A72+1)/ 



(B22) 



Note that these functions are simply related to the general form of (|2.20| ) . 

Then we proceed counting only linearly independent states which do not contain sym- 
metric polynomial factors. For low angular momenta we get these numbers 



AM 


1 


2 


3 


4 


5 


6 


dim 








1 


1 


2 


2 
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where each state in the table is S z — element of a triplet, S = 1. The total number of 
low-lying fermion edge states in the untwisted sector of the HR state is then as given in the 
table in §IIC. 

Finally we note that the sets of linearly independent functions obtained here for each 
AM can be rearranged into the general form derived in Appendix A. 
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